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ABSTRACT 


Orientation and growth of industry towards batch production has spurred research into 
developing efficient means for solving its core design and operational decision 
problems. The result is the development of Cellular Manufacturing System (CMS), 
which evokes the appUcation of Group technology. Cellular manufecturing (CM) 
involves processing set of similar parts on a dedicated cluster of machines. One of the 
problems in the design of CMS is the identification of part families and machine cells, 
which leads to simultaneous production of a family of related part in a flexible 
manufacturing system (FMS). This partitioning of the manufacturing system helps 
decompose the FMS planning problem into more manageable and more autonomous 
sub-units. 

The aim of the present work is to develop solution methodology for finding part 
famihes and machine cells for some of the commotdy occurring grouping scenario in 
CMS. The emphasis is particidarly on generalized grouping problem where each part 
can have alternative process plans and each operation m these plans can be performed 
on alternative machines. Literature turns up mto various techniques of solving this 
problem. The problem turns out to be NP-conplete, hence a number of heuristic and 
hybrid approaches besides mathematical programming models have been developed to 
obtain optimal or near optimal solutions. In this work, an integer programming model 
named Mutually Adjacent Association Model (MAAM) has been proposed. The model 
is shown to be equivalent to the unit capacitated circulation network flow (UCCNF) 
model. The objective is to maximize (minimize) the measure of similarity (dissimilarity) 
between the process plans of parts and then to group them mto families such that all 



the parts referred by the plans selected, can be processed in their correq)onding 
machine cells. The model does not require any knowledge about the number of part 
families to be formed. The solution to the problem is obtained in terms of one main 
loop and several sub-loops identified by using CPLEX. The concept of applying 
network flow in formulating the problem as a mathematical model is relatively a new 
approach. A simple hemistic is proposed for getting the final solution of part famibes 
with their corresponding machine cells. It appears that the proposed model (MAAM) 
has excellent potential to provide computationally optimal solution conq)ared to the p- 
median model (PMM) by Kusiak[1987] and Agarwal’s modified model. 

Some measures of similarity and dissimilarity between the pairs of process plans have 
been proposed which are able to capture the similarity or dissimilarity between the 
plans with respect to processing time of operations and their sequence of occurrence. 

Besides these, improvements in PMM have also been modeled which have the 
flexibOity of either accepting the number of famibes in advance, or escaping Jfrom it. In 
the former case, termed as inq)osed grouping model (IMPGRP), one needs to inquire 
about the maximum number of parts to be kept in a family so that the conq)utational 
time is the least; while in the latter case, termed as natural grouping model (NATGRP), 
the in:q)roved model furnishes with the number of part famibes by itself and one need 
not bother about the maximum number of parts to be kept in a family for minimum 
computational time to result. To obtain the final optimal solution, the hemistic 
proposed is appbed. Computational results shows that NATGRP is much faster than p- 
median model (PMM) while IMPGRP is weaker in this sense. 
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INTRODUCTION 


Group Technology (GT) is the realization that amidst diversities many problems are 
similar, and that by grouping similar problems, a single solution can be found to a set of 
problems thus saving time and effort. 

In the context of manufacturing, GT is an organizational philosophy that aims to collect 
manufactured parts and components into groups, based on their similarities, to facilitate 
their production and at the same time resulting into the effective use of manufacturing 
resources. This is essentially the forming of groups of parts (or part families) with s imilar 
characteristics from a design or manufacturing point of view to reduce the number of 
unique tasks to be dealt with. The idea behind GT is, thus, to deconq)ose a manufacturing 
system into subsystems. In addition to the simplification of management control through 
the creation of subsystem called cell, GT leads to higher manufecturing efficiency and 
many other hidden benefits. A cell is essentially a group of machines dedicated to a family 
of related components. The cell is self contained and self managed, with performance 
control entrusted to the foreman v^dio wields authority to change loading on machines, 
assigns operations to machines or alter the sequence of operations that he deem fit for 
proper load balancing within cell so long as such decisions do not contradict the overall 
strategies by top management. 

Group Technology has probably had greatest impact on increasing manufacturing 
productivity than any other manufacturing concepts. It is recognized as an essential 
element of the foimdation for the successful development and hrplementation of 



automated manufacturing system. This is due to the contributions the cellular 
manufacturing concept has made to other manufacturing technologies and approaches 
such as Computer Aided Process Planning (CAPP), Computer Aided Manufacturing 
(CAM) and thereby to roboticised and flexible manufacturing systems. Recognizing the 
contribution that GT can make to the implementation of CoiDputer Integrated 
Manufacturing System, Merchant says, "the appropriate initial step, to lay a sound 
foimdation for the gradual evolution of a factory to full computer control, is to institute 
group technology cellular organisation" (Choobineh, 1988). The apphcation of GT does 
not depend upon the degree of automation, and may be applied in a totally automated or 
even in a manual production system (Choobineh, 1988). 

1.1 DESIRABILITY OF GROUP TECHNOLOGY 

In the conq)etitive market, an industry involved with small lots and a large variety of parts 
covets to be more and more productive and at the same time more flexible to respond to 
the call of new and diverse products in the least time. As we know, the ideal production 
liue for any product is mass production type wiiich has a flow line arrangement J&cihtating 
in smooth flow production and low work-in-process (WIP) inventory. Batch and job 
production, however, are generally forced to adopt functional layout with a large number 
of general purpose machines to suit the economies of medium volume and low volume 
production. Consequently, such systems are afilicted by problems of long queues in fi'ont 
of facilities, unnecessary setup changes between parts, large amount of material handling 
and long manufecturing lead times. 

These disadvantages and uneconomical productivity are the main concerns in the 
emerging business world. And the response of the newer manufacturing technologies, 



both of the machines and of the tooling, to mitigate them remains inconiplete without an 
at par improvement in organization philosophies because technological advances such as 
improved tools with higher cutting speed, longer tool life alone will not eliminate 
inefficiency since only five percent of the total time is spent in actual processing in 
producing a part. Moreover, technological improvements, are expensive. The 
organizational inprovement is, however, more affordable. One such improvement in the 
organization of production is the appKcation of group technology. 

Group technology with the advantages (Kusiak, 1990) of reduced production lead time 
(20-88%), reduced work-in-process (up to 88%), reduced labour (15-25%), reduced 
tooling (20-30%), reduced rework and scrap materials (15-75%), reduced set-up time 
(20-60%), reduced order time defivery (13-136%), inproved human relation and reduced 
paper work addresses the difficulties &ced by today’s manufacturing industry. Moreover 
the need for flexibility and just in time processing capabilities is forcing many discrete part 
manufacturers to consider reorganizing their fecilities. 

1.2 SCOPE OF PRESENT WORK 
(a) Mathematical Formulation 

Group technology problem has been of tremendous interest for practitioners as well as for 
academicians during the past three decades. Numerous analytical or heuristic methods 
have been developed to solve the problem of forming machines into groups and 
components into associated families or vice versa. But there has been httle effort in 
developing methodologies applying network flow approach, an efficient mathematical 
took It has been pointed out by many researcher that grouping is non-polynomially (NP) 
hard problem. A pioneer work, perhaps the first of its kind using this approach has been 



developed by Lee and Gracia-Diaz £1993] wiiich formulates the sirtq)le grouping problem 
as unit capacitated circulation network with the distance measures between any two 
machines as the per-imit flow cost of capacitated network. Another graph theory 
approach based on network partitioning has been applied by Rajagopalan and Batra 
[1975], Vohra et. al.[1990] and Wu and Salvendy [1993], All these studies using graph 
theory assume that each part has only one process plan and each operation of the process 
plan can be performed by only one machine. Generalized grouping concept, where a part 
can have more than one process plan due to the presence of alternative machine, was 
introduced by Kusiak [1987]. The group formation problem is to select a speciBc process 
plan for each part and group them into part families. He proposed a 0-1 hnear integer 
programming solution known as p-median model (PMM) for this problem which assumes 
that number of groups is known in advance. But, this is computationally expensive even 
for a meditim size problems (Wei and Gaither, 1990). Also, it is executed for a number of 
trials with different number of families to reach at the optimal solution. A more general 
problem environment is where each part has more than one process plan and each 
operation in a process plan can be performed on more than one machine. Rajmani 
et.al.[1990] proposed general linear and nonlinear integer programming formulations for 
such problem scenario which are NP hard in strong saise (Logendran et. al, 1994). 
Kusiak and Cho (1992) suggested an algorithm to solve part family formation in 
generalized grouping by solving for the cliques in a graph of nodes of process plans. This 
method, however, can solve the grouping problem only when the digoint families exist 
which, in general, is not the case. Nagi et. al.[1990] proposed an algorithm which 
considers multiple process plans along with capacity of machines. The objective is to 
minimize the intercell movements. The algorithm solves the problems of process plan 



selection and machine ceU formation iteratively. Won and Kim [1994] proposed an 
assignment model to solve cell formation with multiple process plans. The algorithm is an 
iterative process using a heuristic procedure. Motivated by the work of Lee and Gracia- 
Diaz [1993], a new integer programming formulation named Mutually Adjacent 
Association Model (MAAM) using the apphed concept of capacitated circulation network 
is being proposed to solve the generalized grouping problem. The formulation does not 
need the niunber of part families to be prespecified. However, as the formulation results 
into many families which are not disjoint, a heuristic procedure is proposed to merge two 
or more than two of them. The proposed heiuistic also involves machine cell formation 
with the objective of minimizing the intercell movements. This is not an iterative process. 
Once the problem is solved it will give the optimal solution in one go. The proposed 
formulation along with formulation of generalized p-median method (Kusiak, 1987) and 
modified Kusiak’s model (Agarwal, 1990) have been coded and applied to solve some 
examples from literature. Since the order of complexity (in terms of total number of 
enumerations of forming femilies) of the proposed formulation is very less (of quadratic 
order) as copperas to the p-median and Agarwal’s modified model (of exponential order) 
it has been found to be faster than p-median formulation and gives better grouping in 
terms of exceptional elements. It has also been found that Agarwal’ s modified model is 
not able to solve even a middle size problem of 20-25 plans. 

Improvements in the Kusiak’s model (PMM) have also been proposed and tested. While 
in the Kusiak’s model it is possible that a single part can make a family, thereby violating 
the GT concept itseh^ these improvements agree with the GT concept. These 
iiq)rovements are such that the solver have the flexibility in either prespecifying the 
number of part families or he may escape from this. They have been termed as imposed 



grouping model (IMPGRP) and natural grouping model (NATGRP) respectively. 
NATGRP gives the number of families as its result. However, since many families 
resulting are not disjoint or have common machines as their requirement, the proposed 
heuristic is applied to reach at the optimal solution. A con 5 )arison of these inqDrovements 
with the Kusiak model have also been carried out which diows that NATGRP is 
computationally very efiBcient as compared to Kusiak’s model and IMPGRP. 

(b) Measures of Similarity and Distance 

Similarity (distance) coefl&cient (SC), a numerical measure of commonahty (dissimilarity) 
between macliines or parts, has been used in most of the methods known so far to solve 
the clustering problem Shafer et. AL [1990] presented a survey on this and developed a 
new similarity coefficient. Mosier [1989] ejq)erimented with some well known SCs 
designed for applyiag hierarchical cluster analysis to the GT machine ceU formation. Very 
few SCs have been developed to represent the operation sequence similarity. We propose 
a dissimilarity coefficient, named as Operation Sequence and Machine Utilization 
Dissimilarity Coefficient (OSMUDC) between the process plans taking operation 
sequences and utilization penalty of machines into account. Some examples from the 
literature of grouping are solved with this dissimilarity coefficients. Results are con^ared 
with that of resulting from the application of some existing SCs on operation sequence. 
The coefficient simultaneously offers the most smooth machine layout in the cell, as it 
minimizes the material flows by considering the operation sequence. This coefficient also 
provides the optimum linear machine layout in a machine cell In addition to this, two 
other measures related to process time have been devised. One called the Multiplicative 
Weighted Similarity Coefficient (MuWSC) ensures the minimization of intercell 
movements of parts while the other called Time Movement Distance Coefficient (TMDC) 



ensures the mmimization of time of the intercell movements. Another measure weighted 
part and time similarity coefficient (WPTSC) by makmg a tradeoff between the two is also 
proposed. 


1.3 ORGANIZATION OF THE THESIS 

Chapter 1 furnishes introduction of the subject of group technology. Desirabihty of the 
group technology has also been discussed. 

Chapter 2 presents a brief overview of the major efforts in group technology. Moreover, 
discrepancies and observations are also presented. 

Chapter 3 describes some selected generalized groupiag techniques. Their comparison has 
also been performed through a table. 

Chapter 4 presents the proposed mathematical model (MAAM) and the improvements in 
the existing Kusiak’s model (PMM). Comparison of the proposed model and the 
improved models with Kusiak’s model has been discussed on the basis of a number of 

V 

criterion. 

Chapter 5 devotes to the similarity and distance measures developed. Also, some existmg 
measures parallel to them have been discussed. 

Chapter 6 explains various efficiency measures to compare the results of various grouping 
techniques. 

Chapter 7 implements and experiments with these models. It discusses and conpares then- 
results for some exanple problems. The graphs, charts, and tables used to analyze the 
solution are presented in the chapter. 

Chapter 8 concludes and brings to view some further scope of research in the related 


field. 



2 

PREVIOUS WORKS : A SYNOPSIS 

In this chapter a brief description of some of the major efforts in the field of group 
technology is presented. Since application of group technology aims at dividing the 
production system into subsystems, most of the research is focused on dividing the 
machines and components into groups. There have been numerous application of 
mathematical programming, heuristics, neural networks, fuz 2 y mathematics, expert 
system etc. \vith diOferent kinds of inputs to suit a manufacturing 
environment. 

A combination of inputs or a single input from the list given below is usually taken 
for the purpose of solving the GT problems. 

■ Machine-part 0-1 matrix 

■ Machine-part routing sequence 

■ Process time 

■ Production volume 

■ Alternate part routing 

■ Machine capacity 

■ Costs (Handling, operation, overtime) 

The organization of description of major efforts in the field of GT is based on 
some classification scheme scheme as described below: 

1. Rules of Thumb 

2. Classification and Coding System 

3. Flow Analysis 

Approaches using Similarity/Efissimilarity Coefficients 


4. 



5. Numerical Taxonomy 

6. Mathematical Programming Approach 

7. Heuristic Approach 

8. Hybrid Approach 

9. Network Approach 

1 0. Fuzzy Mathematical Approach 

1 1 . Artificial Intelligence System 

2.1 Rules of Thumb 

Based on the production engineer's local knowledge of the parts and processes, 
some femihes of the components may be self evident. Families consists of 
components defined byname or functions e.g. punrq) casing, shafts and sleeves. 

2.2 Classification and Coding System 

In this system, the first stq) is to classify and code all the con: 5 )onents using a 
stable system where, for every manufacturing and design attributes of the 
components a unique code is reserved. In the production department code wise 
route cards are prepared for all the corc 5 )onents and con:q)onents having similar 
codes are identified to form con 5 )onent families and thereby to form associated 
machine cells. In the design office the code will be useful in retrieving earher 
designs of similar form. The problem of design retrieval is tackled by allocating a 
code number through a classification system which covers the inq)ortant design 
features of the components. 

A niunber of commercial (for e.g. BRISCH, CODE, MICLASS) as well as non- 
proprietary (for e.g. OPTTZ, KC - 1) coding scheme are cuiraitly in practice. 


2.3 Flow Analysis 

The aim of flow analysis is that of finding the families of coniponents and 
associated groups of machines for group layout, by a progressive analysis of the 
information contained in route cards. 

(a) Production flow analysis 

Burbidge [1963] deserves the credit for his pioneering work on group technology 
who developed production flow analysis (PFA). It involves a number of stages of 
analysis which may be briefly described as follows: 

STAGE 1 involves classifying the machines by a number according to type, on the 
basis of operations that can be performed. Machines capable of performing similar 
operations are usually classified with the same type number and the specific needs 
of components for particular machines within the t 3 pe are considered wlien 
allocating these machines to the groups which need them Machines required for 
min or and ancillary operations are excluded fi'om the analysis. Such ancillary 
machines are appropriately assigned o the required machine-conqionent groups 
once these have been detennined firom the analysis. 

STAGE 2 carries the extensive checking of the parts hst and production route card 
information to identify and ensure correctness of the essential information for the 
analysis; namely, for each coirponent, the operations to be undertaken and the 
machine necessary to perform each of these operations. 

STAGE 3 is termed factory flow analysis v^ch, as the title inq)lies, involves a 
macro examination of the flow of conq)onents through the machines, which, in 
turn, allows the problem to be decomposed into a number of major machine- 
component groups. 

STAGE 4 is to undertake the final group analysis on each of these major groups 
separately, rather than to atteicpt it on the undivided overall aggregation of 
machines and con:q)onents. In feet, this has been the focus of attention of several 
researcher \\dio take the necessary data’s available as a result of the earlier stages 
of analysis. These data’s can be of the type of input(s) listed above. 



(b) Con] 5 )onent flow analysis 

Component flow analysis (CFA) was developed by El-Essawy [1971] wdio claimed 
it to be dififerent from PEA, but the similarity of the two approaches is apparent. 
CFA employs three stages of analysis. 

STAGE 1 objective is to consider the total component mix of the company and to 
identify and sort components into categories according to their manufacturing 
requirements. The conqponents are sorted in order of machine requirements and 
sorted list is printed in two ways, firstly in the order of the number of machines 
required and secondly in the order of the smallest machine number involved. 

STAGE 2 aims to obtain groupings of machines using the lists of sorted 
components and taking into accoimt various local constramts. Rough groups are 
formed by using the combinations with the highest number of machmes as the 
cores, to which other machmes and components are successively added. 

STAGE 3 involves a detailed analysis of the loading and flow pattern of the cells 
with appropriate adjustments to ensure that an acceptable design is achieved. 

CFA does differ from PFA in some reg)ects. The latter first partition the problem, 
whereas the former does not. The manner in which the cells are built up is also 
different in the two methods. 

2.4 Approaches using Similarity/Dissimilarity Coefficients 
The similarity coefficient approach was first suggested by McAuley[1972]. This 
class of methods calculate a nmnerical measure which describes the similarity 
between parts or machines. A survey of similarity and distance measures applicable 
to cellular manufecturing has been presented by Shafer and Rogers [1993]. These 
approaches are divided into two two subclasses hierarchical clustering and non 
hierarchical clustering. 

(a) Hierarchical clustering 

The basis of this method is to measure the 'similarity' between each pair of 
machmes (parts) and then to group the machines (parts) into families based on their 


similarity measurements. In most cases, the similaiity measurement used is the 
coefiBcients of Jaccard which is defined for any pair of machines as ‘the number of 
components which visit both machines, divided by the number of con:q)onents 
which visit at least one of the machines'. The con^onents having similar 
coefficients are grouped together to form a family 

McAuley [1972] used Single Linkage Cluster Analysis (SLCA) based on this 
similarity coefficient. This method first clusters together those machines mutually 
related with the highest possible similarity coefficient, then it successively lowers 
the level of admission by steps of predetermined equal magnitude. 

Witte [1980] categorized the roachines as primary, secondary and tertiary based on 
the interdependence of the three classes of machine types and proposed primary, 
secondary and tertiary similarity coefficients. 

Selection of threshold value of SC in SLCA, Witte's three SCs and Rajgopalan's 
graph theoretic [1975] is arbitrary and it requires a certain amount of judgment. As 
shown by de Beer and de Witte [1978], the Jaccard coefficient fails vriien one of 
the machines processes a large nmnber of parts which are not common to the 
machine pairs. 

(b) Non hierarchical clustering 

These methods first select seed machines and then cluster all the machines around 
these machine. Similarly parts are clustered around seed parts. These seed are 
selected firom some mathematical programming approach. ZODIAC 
(Chandrasekharan and Rajagopalan, 1987) calcxilates the seed fi'om the bipartite 
graph formulation, while GRAFICS (Srinivasan and Narendaran, 1992) obtains the 
initial seed firom the assignment method. 

2.5 Numerical Taxonomy 

Numerical taxonomy is the widely used technique for cell formation in group 
technology. These methods re-arrange the rows and columns of the machine 


component incidence matrix such that the I's are brought together. Each block of 
I's constitute a component family and a machine cell. 

McCormick [1972] proposed Bond Energy Algori thm (BEA) which reorders the 
rows and columns of the matrix for the purpose of moving numerically larger 
matrix elements together. A measure of closeness between two corrponent vectors 
or machine vectors termed bond energy is defined. Method seeks to determine 
permutation of rows and columns in which the sum of products of adjacent 
elements is maximized. This is a restricted form of quadratic assignment problem. 
Graham [1976] proposed linear weighting algorithm where the wei^ts are 
increased linearly. The ith row is a given a weighting of m - i + 1 where m is the 
total number of rows, and the priority ranking value is determined as the mean of 
the weightings of the non zero entries. Ranking values calculated this way can be 
found and sorted very quickly and the requirement of a very large integer 
representation does not arise. The major disadvantages of this linear weighting 
algorithm are the complicated and very confusing patterns of the intermediate 
results as well as the difficulty in predicting the behaviour of llie procedure. 

King [1980] proposed a Rank Order Clustering (ROC) technique based on 
numerical taxonomy. This method interprets each entry in the matrix as a binary 
word and confutes decimal equivalent of each row and column. Rows and 
columns are sorted in decreasing order of the binary word values, until algorithm 
converges. 

King and Nakomchai [1982] observed that reading entries as binary word 
restricted the size of the problem to less than 47 machines and 47 components 
because the largest integer rq)resentation in any conputer was (2^^ - 1). They 
proposed ROC2 which is conputationally efficient. 

Chan and Milner [1982] proposed direct clustering \\iiich is a poor version of ROC 
algorithm, except that it eliminates the sensitivity of the latter to the configuration 


of the initial matrix. This algorithm has the same handicaps as the ROC algorithm 
and uses a limited binary comparison procedure for ranlcing. 

Kusiak [1985] discusses the conqtutational complexities of BEA and ROCs. 
Chandrdiekharan and Rajagopalan [1986] observed that ROC does not provide 
block diagonal structime even in the case of structured matrices where such a 
possibility exists. They found that ROC is dependent on the initial disposition of 
the incidence matrix and that it does not handle exceptional elements too well. 
They suggested MODROC. 

Khator and Irani [1987] proposed a heuristic procedure, the Occupancy Value 
Method (OVA), for identifying clusters in a machine component matrix created 
from route card data. A unique feature of this method is that it progressively 
develops block diagonalization starting from the northwest comer of the matrix. 
Kusiak and Chow [1987] proposed cluster identification algorithm (CIA) which 
identifies the clusters in a binary machine-part incidence matrix provided they exist. 
He proposed three algorithms based on cluster identification to solve GT problem 
with bottleneck conrponents and machines. The first algorithm solves an 
rmconstrained GT problem. The second is a heuristic which considers a constraint 
restricting the number of machines in each cell and identifies the bottlenecks. The 
third named the cost analysis algorithm allows one to solve the group technology 
problem with known subcontracting costs. The overlapping parts are removed in 
such a way that the total sxim of subcontracting costs is minimized. 

Boe and Chang [1991] identified several deficiencies of clustering and array based 
sorting algorithms and suggested a clustering algorithm termed as close nei^bour 
algorithm. In the first stage clustering of machines is performed through a 0-1 
quadratic programming model maximizing the closeness of machines and in the 
second stage clustering of parts is done. The algorithm converges in one pass in 
most of the problems. 



2.6 Mathematical Programming Approach 

Using mathematical approaches like linear programming, 0-1 integer programming, 
dynamic programming, goal programming etc., the GT cell formation problem 
may be given concise mathematical formulation and a solution can be found for a 
given constraint set. 

Gunasingh and Lashkaii [1989] proposed two 0-1 integer programming 
formulations for the GT problem. The first formulation groups machines based on 
the compatibihty of processing components, while the second formulation groups 
machines m order to minimize the cost of machine allocation and the cost of 
intercellular movement. 

Gunasingh [1991] presented a non-linear 0-1 integer programming formulation to 
simultaneously group machines and parts in cellular manufacturing systems based 
on tooling requirements of the parts., tools available on the machines and the 
processing times. The formulation takes into account the limitations on the number 
of parts and machines in a group and the number of machine types available. 

Boctor [1991] presented 0-1 integer programming based formulation for cell 
formation with the objective of minimizing the exceptional elements. This model 
allows the designer to control the cell sizes. 

Kusiak [1987] developed a integer programming problem. He considered presence 
of alternative process plans and suggested a 0-1 integer program based p-median 
model (PMM) to determine the optimal process plan and assignment of parts to 
machine cells. Incorporating process plans will result in an improved quahty of part 
femilies and machine cells. 

Song and Hitomi [1992] considered the part routing, part volume and cell size. 
They formulated the problem in two steps : first, a QAP formulation to minimize 
the number of parts within a cell and the second, a graph formulation for 
minimization of intercell part flows. Global optimum is obtained using branch and 
bound method. 



RajTnani [1992] developed an mixed integer programming model to minimize the 
saun of total discounted cost and set up cost. He considered dependent set up time, 
cost for switching the product, cost and capacity of machine. 

Rajmani [1994] developed three integer models for the same problem. First model 
assigns machines to parts using any clustering algorithm. Second model assigns 
machines to known part families. Part families are obtained by part attributes. Third 
model obtains part families and machine groups simultaneously considering 
demand time and resource constraint. He considered alternate machines and 
alternate routings for part, and cost and capacity of machine. 

Wei and Gaither [1990] developed a 0-1 integer programming model to assign 
machines and parts to cells simultaneously. This model can solve larger cell 
formation problems than earlier approaches because it uses far fewer decision 
variable than earlier p-median formulation. The objective of this model is to 
minimize the cost of manufacturing exceptional part types outside the cell subject 
to machine capacity constraints. 

Kumar [1986] formulated a joint grouping and loading problem as a multistage 
multi-objective modeL The min-max approach to multiobjective optimization is 
used to obtain a compromise solution. 

Shtub [1989] showed that the generalized cell formation problem due to Kusiak 
[1987a] is equivalent to the generalized assignment problem (GAP). 

Agarwal [1990] investigated p-median formulation of Kusiak [1987] and GAP of 
Shtub [1989] and improved formulations requiring lesser number of variables and 
constraints. 

Choobineh [1988] proposed a two-stage procedure to design a cellular 
manufacturing systenr In the first stage, part family is determined by similarity 
coefficient based algorithm. In the second stage an integer programming 
formulation is used to form machine cells with an objective of minimizing 
production cost and the cost of acquiring and maintaining machine tools. 



Srinivasan [1990] proposed an assignment model for GT. The similarity coeflBcient 
matrix for machines is solved as an assignment problem from which initial cells are 
identified. Coitq)onent families are similarly identified and they are merged using a 
set of rules. 

Dahel and Smith [1993] proposed a multiojective model for grouping parts and 
machines simultaneously into prespecified number of cells. The objective of the 
model is to minimize inerceU moves as well as to maximize the processing 
flexibility by maximizing the number of machines in each cell subject to machine 
capacity and cell size constraints. 

2.7 Heuristic Approach 

Since even the relaxed version of the grouping problem is NP-complete, it is 
unlikely to ever get the optimal solution to the original problem efficiently. 
Therefore, the cell formation procedures found in the literature are predominantly 
heuristic which can find good solutions in a reasonable amount of con^utational 
time. The heuristic can be problem specific or general purpose. 

Lenstra [1972] pointed out that the chisterittg problem is equivalent to two 
travelling salesman problems. The computational conq)lexity of the two problems 
is non-polynomially-con^lete (NP) for large problem 

Shtub [1989] showed that the simple cell formation problem is equivalent to the 
generalized assignment problem (GAP) >\bich is NP-hard as the partitioning 
problem (NP- complete) is reducible to GAP (Fisher et.al., 1986). 

BaUakur et. aL [1987] have proposed a within-cell utilization heuristic for 
designing cellular manufacturing systems. The distinguishing feature of this 
heuristic is its consideration of several practical criteria such as within-cell 
utilization, work load fractions, maximum number of machines that are assigned to 
a cell, and the percentage of operations of a parts completed within a cell 
Vannelli et. aL [1986] discussed a method of identifying the minimal number of 
bottleneck cells vsiiich, when dealt with through either diq)lication of machines or 



subcontracting of parts, will result in perfect part machine grouping with no 
overlap. The polynomially bounded algorithm used in the analysis are oriented 
towards finding minimal cut-nodes in either partition of the bipartite part-machine 
graph. 

Askin et. el. [1991] proposed a hamiltonian path approach to reordering the part 
machine incidence matrix for ceU formation. Similarity coefficients have been used 
to form a distance measure for machines. Using the distance matrix, a heuristic 
solution to the associated travelling salesman problem (TSP) was obtained. Finally 
using the heuristic tour determined for TSP a heuristic solution to the associated 
hamiltonian path problem was obtained. 

Harahalakis et. al. [1990] proposed a simple twofold heuristic algorithm capable of 
minimi ring interceU material movement by considering the sequence of operation. 
The fibrst stage is a bottom-up aggregation procedure to minimize the ‘normalized 
inter ceU traffic’defined as the ratio of the inter ceU movement between any two 
ceUs imder consideration and the number of noachines included in these cells. The 
second stage is aimed at mq)roving the solution obtaiaed firom the first stage, in 
which the significance of a machine to a ceU is validated. 

Logendran [1990] proposed a model which incorporates both mtraceU and interceU 
moves. Logendran [1991] considered sequence of operations iato account hi the 
design of ceUular manufacturing system. This model mcludes the mqiact of ceUs hi 
evaluating interceU moves. The model also aims at maintaining a targeted 
utilization of work station. 

Balasubramanium [1993] described an algorithm which helps the formation of 
machine groups. The algorithm first considers several possible ceUs. It thai 
develops a similarity matrix to reflect the conqionent’s routeing simUarity with that 
of the arrangement of the facilities in the ceU. He developed similarity coefficient 
considering production volume and operation simUarity. The machine grouping is 
done heuiisticaUy by considering the cost of handling (using SC) and the cost of 



idle time of machine (using machine loads) and overtime sunultaneously. The same 
machine type is allowed m more than one cell if necessary. 

Venugopal and Narendran [1992] presented an algorithm based on simulated 
annealing to solve the machine-part grouping problem The algorithm partitions the 
machine-part incidence matrix in such a way that the total intracell load variation is 
minimal. 

Genetic algorithms (GA) a general adaptive search method which mimics the 
process of natural evolution were developed by Holland [1975], Since then GAs 
have been applied to diverse areas of optimization. Its success in combinatorial 
optimization attracted the researchers of GT to apply GA for solving grouping 
problenns. Venugopal and Narendran [1992] seem to be the first who applied GA 
for solving grouping problems. They formulated the grouping problem as a 
bi-criterion mathematical program and solved it using GA. Computational results 
show that the proposed algorithm can be successfully applied to large scale 
grouping problems. 

2.8 Hybrid Approach 

The survey tiU now has shown that certain ceU formation methods employs both 
heuristic and mathematical procedures in different stages. Such methods are called 
hybrid approaches. 

Steudel and Ballakur [1987] proposed a two-stage heuristic in which the first stage 
employs a dynamic programming procedure to generate an optimum chain of 
machines in which the sum of bonds between machines in the chain is maximized. 
The second stage is a heuristic, that partitions the chain to form machine cells 
subject to cell size restrictions. 

Nagi et. al. [1990] addresses two problems in the presence of mult^le routeing 
process plans for a part : (a) process plan selection which has been solved by 
formulating it as a LP problem, (b) cell formation, is solved throu^ a heuristic. 



The method also considers the sequence of operations and production data such as 
production volume, processing time and machine capacity. 

Choobmeh [1988] employed a two stage procediure. First ertq)loys a similarity 
measure based on part operations and operation sequences and attempts to identify 
part families using a heuristic procedure. The second stage formulates as an integer 
programming which specifies the type and number of machines in each cell and the 
assignment of part families to the cells. 

Askin and Chhi [1990] incorporated costs of inventory, machine depreciation, 
machine setups and material handling costs into a mathematical programming 
formulation. A heuristic graph partitioning procediue is then applied to 
solve it. 

2.9 Network Approach 

The network approach to GT problem can be divided into two categories. 

2.9.1 Network partitioning approach 

Machines and parts are treated as nodes and machining time or SC is put on the 
arcs connectmg them. These models aim at obtaining disconnected subgraphs for 
identifying cells and part family. 

Rajagopalan and Batra [1975] suggested a graph partitioning approach vdiich used 
cliques of machine-graph as a means of grouping machines. He used the Jaccard 
SCs on the arc. 

Vanneli and Kumar [1986, 1987] have developed algorithm to identify minimal 
bottleneck machines to be duplicated and components to be subcontracted to 
obtaia non-intersecting part femilies and cells for a given matrix. 

Vohra et. al [1990] preesented a non-heuiistic network approach to form machine 
cells with minimum iatercellular interactions. The machme part matrix containing 
machme times is represented as a network which is subsequently partitioned by 
using a modified Gomory-Hu algorithm to find a minimum cut between nodes of 
machines. 



Wu and Salvendy [1993] proposed a modified version of Vohra's work and for 
larger problems it gives a good solution but not the optimal solution. 
Graph-deconqjosition is a difficult combinatorial task and the above mentioned 
methods do not guarantee optimahty. 

2.9.2 Network flow approach 

Lee and Diaz [1993] have described a simple groupmg problem as a maximum 
flow network problem. Each machme is represented by two nodes and each arc has 
three values (U,6ij,L) assigned where 5ij is the SC between machine i and machine 
j and U and L are the upper and lower limit on the flow respectively. 

There has been no work foxmd to solve generalized grouping problem using this 
approach. 

2.10 Fuzzy Mathematical Approach 

Most of the methodologies for cell formation discussed earlier assume that the 
information about processing costs, processing times, part demand, etc. is precise 
(deterministic) but in reality it is not. Fuzzy clustering provides a solution to such 
problems. Relatively few papers such as Xu and Wang [1989] and Chu and Hayya 
[1991] have addressed the issues of vagueness in cell formation problems. 

2.11 Artificial Intelligence (AI) Based Techniques 

Recent development in AI techniques for solving decision making problem in 
various fields have influenced manufacturing and group technology, too. A brief 
review of the AI based approaches for single grouping is given below. 

2.11.1 Expert system and knowledge based system 
Kusiak [1988] first proposed EXGT-S, the tandem architecture of knowledge 
based system for simple grouping problem. It has two parts, namely expert system 
and heuristic clustering algorithms. Each partial solution generated by heuristic 
algorithm is evaluated by the expert system which modifies search direction of 
algorithm. The approach presented allows one to take advantage of the user's 



production expertise, and can be used to determine machine cells and part Emilies. 
However, the EXGT-S is unable to detect the machine chaining problem. Chow 
and Hawaledika [1993] presented another knowledge based system which can be 
used for identifying machine cells and part famili es and helps in rectifying some of 
the problems commonly encountered in a EXGT-S system 
2.11.2 Neural network based methodologies 

The nemal network is a new algorithmic approach which has been the subject of 
intensive study of mathematician, statisticians, physicists, engineers and conqjuter 
scientists. Its parallel processing capability, learning ability, and less storage space 
requirements have made it popular for solving large scale real life problems. Kusiak 
[1991] proposed adaptive resonance theory based neural network developed by 
Carpenter and Grosberg to solve the GT problems. 

Karpathi and Suresh [1992] used Carpenter and Grosberg neural network. Since 
this clustering method utilizes binary-valued inputs and it can be trained without 
supervision. Experiment on larger data sets (10000 x 100) revealed that the 
method results m the identification of clusters with fast execution times. 

Chu [1993] applied a neural network approach based upon a competitive learning 
paradigm for designing part famihes and machine cells. The procedure can. 
effectively obtain optimal cluster result. 

2.12 Observations 

The most important point to be noted is that optimum result of GT problem may or 
may not be possible to get because of its NP-complete nature. Basically all the 
research papers have moulded the original GT problem according to some 
constrained environment. Original GT problem is not just the grouping on the 
basis of machine part matrix, rather it should get integrated with material 
requirement planning, just-in-time, market demand structure, assembly operations 
and many more which influences the operation of the production shop. 



It is because of its NP-conqjlete nature that it is not even possible to predict 

whether an optimal solution to a problem exists 

Based on hterature survey presented following observations can be made ; 

• Less efforts have been made to incorporate flexibility and its measure in the 
design of group technology. 

• Network flow approach, an efficient mathematical tool, has not been apphed to 
the generalized groupmg problems. 

• Less work on facihties design (layout) along with cell formation with required 
sequence to facilitate smooth flow production. 

• In most of the approaches production volume and the machnie loads were not 
considered. In some case it was considered only after the cell formation was 
over. 

• In all the SC methods involving production volume and workloads, SC is 
mdependent of production volume. 



SOME TECHNIQUES FOR 
GENERALIZED GROUPING 

This chapter presents the details of some of the techniques which are used to solve 
the GT problem where every part has a set of process plans. However, all of them 
can not be equally optimal m respect of various manufacturing yardsticks like total 
time taken by the plan, number of types of machine required, skill requirement, and 
of course the most important being quahty control. But with reference to the 
concept of GT, one might even settle to choose a sub-optimal plan in order to 
confine a part q)ecific cell. This may be preferred to reduce queuing times, 
transportation times, costs, and WIP, and to increase productivity. Another reason 
for alternate process plans is the existence of fimctionally similar workcentres i.e. 
alternative machines. In a functional layout, this is not a consideration, because the 
part can be routed to any such available workcentre. For example, a part requires a 
turning operation. Our turning facihty consists of two identical lathes. The part 
could use either lathe when it is routed to the turning facility. On the other hand, in 
a GT environment, manufecturing cells usually consists of functionally dissimilar 
workcentres; the two lathes are hkely to be placed in different cells. In this case, 
we would prefer to route the part to the lathe in its corresponding cell Thus, there 
is a need to identify each workcentre as a specific one. This unique identification of 
functionally similar workcentres lead to the existence of mult^le routeing. 

Thus, the cell formation problem in the presence of alternate routeing is to choose 
one of them to determine the part families with most similar route of other parts. 
Then subsequent machine allocation to these femilies produces the machine cells. 
We have selected some of the rq)resentative mathematical approaches which are 
described in the following sections. 



3.1 Generalized p-median method (PMM) pCusiak, 1987] 

3.2 Modified p-median model [Agarwal, 1990] 

3.3 Shtub’s Assignment Model [Shtub, 1987] 

3.4 Maximum Clique Algorithm (MCA) [Kusiak and Cho, 1992] 

3.1 Generalized p-median model (PMM) 

Kusiak [1987] introduced the concept of generalized grouping. Based on this 
concept the part-machine matrix is replaced by the machine process plan binary 
matrix, with the basic and alternate process plans. The objective is to form groups 
of the process plans such that the total similarity (distance) defined in some 
suitable manner is maximized (minimized), and only one process plan from the 
several alternatives of a part is selected in the group formation. This model 
requires the desired number of part families in advance. 

Notation and definitions 
Indices 

i, j : process plans 
k : part 
m : machine 
Parameters 

M : the total number of machines 
K: the total number of parts 

p ; the ma ximum number of process families required 
q : the total number of process plans over aU parts 
Fk : set of process plans for part n 

a;^ = 1 if process plan i requires machine machine m 
= 0 otherwise. 

Sij = similarity between process plan i and j 
The elements of the similarity coefficient matrix [Sy] are computed as follows. 



(a) Similarity between two process plans of the same part 

Sij = -00 (3.1) 

(b) Similarity of a process plan with itself 

Sij = 0 (3.2) 

(c) Similarity between two process plans of dififerent parts 

M 

Sij = Z d(auii,ajni)= V i e Ft ; V j g Ft and k = 1, , K 

m=l 

(3.3) 

Decision Variables 

A process family is identified by the index of one of its member plans. The p 
required families, therefore, will assume their indices fi:om the range 1 to q. 

For i, j = 1, , q, the following decision variable is defined: 

Xij = 1 ifprocess plan i belongs to process family) 

= 0 otherwise. 

Formulation 

The objective of the generalized p-median model is to form the process plan 
famili es by maximizing the value of similar ity measure among them 

q q ^ ^ 

Maximize Z Z siixii w 

i=lj=l •' •' 

subject to 

(I) Selection of only one process plan for a part 

S Zxii = 1 (3.5) 

ieFkj=l 

(n) Formation of maximum required number of groups 

q 

Z xjj < p (3.6) 

j=l 

(in) Assignment of processs plans to groups 

xij < xjj i,j = l, ,q (3.7) 


(IV) Bhiaiy decision variables 
jqj = 0 or 1 


(3.8) 



3.2 Modified p-median model [Agarwal, 1990] 

Agarwal [1990] modified the constraint of assigning the process plans to group in 
generalized PMM. He observed that the number of assignments along the diagonal 
of the similarity coeflScient matrix [Sij] is equal to p (see eqn. 3.6). Thus, the 
remaming (K-p) assignments have to be made coTreq)onding to some non-diagonal 
elements of the matrix in a manner such that tbese (K-p) assignments appear in 
only those p columns or p rows that represent the process families and have 
assignments at their cells corresponding to the diagonal of the matrix. These 
observations have been used to sin:q)lify the process plan assignment constraint 
(i.e. eqn. 3.7) of the PMM. The surq)lified constraint in place of equation (3.7) is in 
the form of following mequality 

q 

Z Xij ^ (k -p +l)xjj j=l, ,q (3.9) 

i = 1 

The idea behind the constramt is that out of assignments of K plans (K-p) 
assignments have to be made from the non-diagonal side and since to any process 
family of p process frmilies to be made, at the most (K-p+1) plans can be assigned 
when only one plan (i.e. the index plan) is assigned to other (p-1) families. 

3.3 Shtub’s Assignment Model 

The grouping problem has been formulated as a generalized assignment problem 
(GAP). The structure of GAP as described by Shtub is as follows. 

Notation and definition 

I set of tadcs 

J set of agents 

i e I index set of tasks 

j e J index set of agents 

Cjj the cost incurred if task i is asagned to agent j 

ly > 0 the amount of a resource required by agent j to perform task i 



aj > 0 the minimum amount of the resouxce that may be expended by 


agent j 

bj > 0 the maximum amount of the resoxnce that may he expended by 
agent j 

Xij = 1 if task i is assigned to agent j 
= 0 otherwise. 

Formulation 

Minimize 

S Z cijxij (3.10) 

i elj eJ 

subject to 

aj^ Z lijxij^^bj foralljeJ (3.11) 

i el 

2 xij = l foralliel (3.12) 

j eJ 

The objective function minimizes the cost of assigning tasks to agents. Constraint 
set (3.11) ensures minimum and maximum usage of the resource by each agent, 
and set (3.12) ensures that each task is asagned to exactly one agent. He also 
showed that the GAP formulation is equivalent to generalized PMM by taking Cy = 
-Sij, and Cjj = -M, a negative number vshere M > niax{Sij} . 

3.4 Maximum Clique Algorithm (MCA) 

Kusiak and Cho [1992] proposed this algorithm to solve the part femily 
formation problem in presence of basic and ahemative process plans of parts. This 
procedure uses the similarity coefficient approach. A similarity coefficient between 
process plans is defined and a branch and bound algorithm is proposed to cluster 
the process plans according to their sirnilarity value. 



The algorithm represents the process plan similarity matrix by a transition graph 
and finds the maximum clique of the graph which is considered as a part femily. 
Similarly all the part families are generated by branch and bound procedure. 
Defiboition 

Sij = similarity coejfficient between two process plans i and j of the process 
plan -machine matrix [aim] is given by 

Sij = 1 if aim > ajm or aim < ajm for all machine m (313) 

= 0 otherwise. 

The value of similarity coefficient Sij = 1 indicates that the one process plan is a 
subset of another process plan. So the parts represented by the process plans do 
not require different type of machines and can be put in a single family. 

Algorithm 

The algorithm works on the branch and boxmd procedure. Similarity coefficient 
between process is calculated and stored in a similarity matrix which is used to 
construct a transition graph. Every node of a graph represents a process plan and 
an arc between two nodes suggests a similarity coefficient of 1 between process 
plans. The disconnected nodes have the similarity of 0 value. All the closed tour 
with the maximum number of arcs between nodes known as maximum clique is 
found out and treated as a process plan families. After deleting the parts 
corresponding to these process plans the similarity matrix is revised and the 
procedure is repeated by the branch and bound procedure resulting into digoint 
part families. 

Steps 

STEP 1: Begin with the generalized matrix [ajm] between process plans and 
machines and set the level L = 0. 

STEP 2; Initialization, Transform matrix [aim] to similarity matrix [sij]. 

STEP 3 ; Represent the current similarity matrix with transition graph. 



STEP 4; Determine maximum cliques in the transition graph hy using any available 
heuristic. Assign level number L= L+1 to each chque and branch on a cHque with 
the maximum total sum of node degree. 

STEP 5: Delete jfrom [sij] the rows and columns correq)onding to parts that have 
process plans included in chque selected. 

STEP 6: Exclude from further consideration a chque which produces larger 
number of process famflies than the current solution. [Fathoming] 

STEP 7; Stop when the current solution is acceptable, otherwise backtrack to 
nearest unfathomed node at level L and repeat from Step 3. 

3.5 Discussions 

In some practical cases maximum number of fenhhes (p) in PMM might be the 
result of experience or planning requirement. If the number of process famili es p is 
not known one can start with a smaU nmnber of process famihes p = Po and solve 
the model. If the quahty of process famihes is not satisfactory p is increased by 1 
and the model is solved again tih a satisfactory result is obtained. 

Table 3.1 presents the details related to the number of variables and constraints 
used by the above models. From the table, it can be sera that for aU practical 
grouping problems, the number of constraints in the Shtub GAP model will be less 
as compared to that of PMM. However, the number of the decision variables wiU 
always be more for GAP. The Agarwal’s modified method requires the least 
number of constraints compared to both the models PMM and GAP and involves 
the same number of variables as for PMM. Thus, the Agarwal’s modified method 
would contain for sure, the least number of constramts and decision variables. 



Table 3.1 Number of constraints and variables for different models 



Kusiak’s 

Agarwal’s 

Shtub’sGAP 


PMM 

modified 




PMM 


Number of constraints 

q^ + K + 1 

q + K+1 

4q + K+l 

Number of variables 



q^ + 3q + 2K 


The most computationally conq)lex steps in MCA are stepl and step3. The 
computational time complexity of step 1 is C>(M^q) where q is the total number of 
rows in the incidence matrix i.e. total number of process plans. The conaplexity of 
step 3 is dominated by clique problem which is NP complete. It’s complexity 
depends upon the heuristic used. So this algorithm is the most sioaple algorithm 
available for the part family formation m presence of alternative process plans. 



4 

SOME PROPOSED MATHEMATICAL 
MODELS 

Bearing the problem environment discussed in the previous chapter, we propose a new 
integer programming model for the generalized grouping problem This problem has 
been formulated as an integer programming model (PMM) by Kusiak [1987] and as an 
assignment model (GAP) by Shtub [1987], Shtub [1987] also showed that his 
assignment model is equivalent to the Kusiak’s integer programming model, however 
there are lesser number of constraints hut larger number of variables (Table 3.1). In the 
hterature of generalized grouping Kusiak’s model is supposed to be the most 
representative mathematical model. It has been foxmd during the survey that most of 
the techniques developed till now have been con^ared with PMM for their 
performance evaluation. It is so because it is the first technique in this field and also 
quite efiScient. But for larger problem the performance worsens in respect of time 
spent in getting the solution. It is because of the complexity of the model being of 
exponential order. It also requires the number of process femilies to be specified in 
advance. This is what we call in:q)osed grouping. This will surely affect the process of 
natural grouping. The satisfectory result is achieved by playing with the number of 
process femihes. His assumption that for the process family formation at least one 
process plan is necessary is also not a practical assumption firom the point of view of 
GT concept. 

(1) The model proposed is inspired fi:om the Lee and Diaz [1993] v^o applied 
network flow approach to solve sinqrle grouping problem vAere each part has a 
single plan. He performed the grouping of machines by minimizmg the dissimilarity 
between machines. The proposed model forms the part families in the presence of 
alternate process plans for a part. Of the several process plans of each part only 



one process plan is selected vdiile forming the part families with the objective of 
maximiTing (minimizing) the similarity (distance) among them. It does not require 
the number of part families to be specified in advance, rather it results in natural 
grouping i.e it gives the number of part families and the part families by itself He 
solution methodology used for the proposed model has been discussed with an 
example problem. Interpretation of the resulting solution which forms process 
famihes is made through the loop formation The model is also described |p be 
equivalent to the Unit Capacitated Circulation Network Flow (UCCNF) proMeiR 

(2) An improvement in the Kusiak’s model has also been proposed. The modified 
formulations have the flexibility of either not ^ecifying the number of process 
femilies in advance or can be specified if so desired by adding the corresponding 
constraint in the formulation Besides it has certain other advantages over both 
the Kusiak’s model and its modified version by Agarwal. A comparison of 
Kusiak’s model, Agarwal’s modified model, the proposed model and the 
improved Kusiak’s models has bear made on the basis of the number of 
constraints, number of variables, order of the complexity, the way of grouping 
and the agreement with GT concepts. 

(3) A sunple hemistic is proposed for the merging of two or more process femilies and 
the assignment of machines to the process families resulting into machine cell 
formation 

In the following sections we describe the proposed model and the inqrrovemait in the 

Kusiak’s model 

4.1 PROPOSED MATHEMATICAL MODEL : MUTUALLY 
ADJACENT ASSOCIATION MODEL (MAAM) 

Definitions and Notation : 

The following definitions and notation are used for the model development: 

M : the total number of machine types 
N ; the total number of parts 
q : the total number ofprocess plans over all parts 



L : the total number of process families 
Indices: 

i, j ; process plans 
k : part 

1 : process family 
m,n : machines 
c : machine cell 

Fk ; set of process plans for part k 
I Fk I = the total number of process plans for part k 
Parameters: 

aim = 1 if process plan i requires machine m 
= 0 otherwise 

Xij = 1 ifplani associates with plan j 
= 0 otherwise 

Umi = Z aims the usage factor for machine m in the process femily 1 
iel 

dij = distance measure between process plan i and j 
Sij = similarity measure between process plan i and j 

It may be noted that process families ultimately corresponds to the part families. 

4.1.1 Distance (Similarity) Parameter Computation 
A distance coefficient between process plan i and j is defined as (Lee, 1981) 
dij = distance measure between process plan i and j 
= ZS(ajjn, ajjn) ; i6Fk,j«?Fk m=ltoM 

= 0; ifi=j 

= +oo; when i,j eFk. (4.1) 

where 

S ( aim, ) " 1 > ^ %n ^ ^'m 
= 0 ; otherwise. 



Similarly, 


Sij = similarity measure between process plans i and j 
“ 2 S ( Ejjjj, Ejjjj ); i eFk, m=ltoM 

= 0; when i=j 

= -co; when i,j eFk. (4.2) 

where 

^ ( ^im’ ^jm ) ~ 1 5 if ~ ^jm 

= 0 ; otherwise. 

Decision variables: 

Xij = 1; ifplani associates with plan j 
= 0; otherwise. 

The square matrix of distance D = [dy] or similarity coefBcient S = [Sy] between the 
process plans thus created following the above definitions is used in the proposed 
model. A distance coefficient is an indicator showing the degree of dissimilarity 
between a pair of process plans. The process plan-machine incidence matrix is used for 
this purpose. In this matrix each machine is identified uniquely by a column number 
and each process plan by a row number indexed consecutively over all parts. The 

matrix A of size q xM is represented by A = [aim], i = 1, , q, m “ 1, , M. 

4.1.2 Integer Programming Formulations 

An integer programming formulation for the process family formation in the presence 
of alternate process process plans is presented below. This formulation assumes that 
the number of process families to be formed is not preq)ecified and that each of the 
process families to be formed has at least two process plans of two different parts 
which, in fact, is in agreement with the group technology concept. 

Objective function 

The objective fimction is to minimize (maximize) the total dissimilarity (similarity) 
between the process plans and can be written as: 



mminwze 


q q 
q q 

or maximize £ Z sij xij (4-3) 

i=lj=l 

Constraints 

1 . Selection of one and only one process plan of a part for association with only one 
of the plans of some other parts, i.e. 

Z Z xij = 1; for all parts k=l, , K (4.4) 

ieFkj^Fk 

2. Selection of a plan not belonging to the part considered (i.e. part k) which 
associates with a plan of the part considered. This plan and the plan selected 
by the constraint 1 (eq. 4.4) belong to different parts other than the part 
considered. 

Z Zxji = 1; for all parts k= 1, , K (4.5) 

j«Fk UFk 

3. part considered the process plans selected by the above two constraints are the 
same so that a Ensuring that for the group formation results. 

Z 2^ xij + Z Zxji = forplani=l, ,q (4.6) 

ieFkjeFk 

It ensures that for a part either a plan will get selected (right hand side = 2), or will 
not get selected (right hand side = 0). 

4. Integrahty of the decision variables 

Xij = 0 or 1; for plans ij =1, ,q (4.7) 

The constraints are ejqplained for example problem given in Table 4. 1 and are diown in 

Appendix A. 



Further, the total number of constraints = K of (4.4) + K of (4.5) + qof(4.6) 

= 2K + q. 

4.1. 2.1 Equivalence of MAAM to the unit capacitated circulation 
network flow (UCCNF) model 

We discuss first (section 4.1.2.1a) the UCCNF model for the machine cell formation 
proposed by Gracia and Diaz [1993] and then we propose the procedure for 
constructing the circulation network for the case of generalized grouping problem 
(section 4.1.2.1b). Gracia and Diaz [1993] solved the problem for the case of simple 
groupmg problem, each part having its basic process plan only. Later we ■will show 
how the proposed integer programming model for generalized groupmg can be viewed 
as an UCCNF model for minimum cost flow problem (MCFP) by analyziag its 
variables, objective fimction and the constraints (section 4.1.2.1c). 

4.1.2.1a UCCNE model for machine cell formation 
(Lee and Diaz , 1993) 

Lee and Diaz [1993] proposed the unit capacitated circulation network flow (UCCNF) 
approach to solve the machine cell formation problem for the case of smq)le grouping, 
each part having their basic plan only. He solved the machine cell formation problem in 
place of part family formation. The argument is that generally the number of machines 
is less than the number of parts. The problem is formulated as the minimum cost flow 
problem (MCFP). The distance parameter between the machines is used as the per unit 
flow cost. 

A capacitated circulation network G = (V, E), where V is the set of nodes and E the 
set of arcs, can be formulated on the basis of the distance parameters between the 
machines. In the network, each machine m is represented by two nodes nia and mt, 
connected by a directed arc (nia, njb). It is assumed that the number of cluster is 
unknown and that each cluster has at least two machines. The procedure followed by 
Lee and Diaz [1993] to construct the network is described below: 

1. Create 2M nodes represaated by ma, mb, m = 1, , M. Arbitrarily, let node 

Ibbe the source node and Lthe sink node of G = (V, E). 



2. Structure the rest of the network according to the following rules. Each arc is 
assigned three values U, L and C to indicate an upper bound on its flow, a 
lower bound on its flow and the per unit flow cost, respectively. 



FIGURE 4.1: A Network Configuration of an M-Machine austering 
Problem 




2. 1 For each node nia , m -i* 1, create an arc directed from ma to tOb with capacity- 
cost triplet defined by [U, L, C] = [1, 1, 0], 

2.2 For each node nib , m =2, , M, create (M-2) arcs directed from mb to ma , 

m = 2, , M , with a capacity-cost triplet [U, L, C] = [1, 0, dn^,] , m n. 

3. For each of the nodes la and lb, create arcs directed from lb to la with 

capacity -cost triplet [U, L, C] = [1, 0, dij, n 1, n =2, , M ; and arcs 

directed from mb (m 9^ 1) to la with capacity-cost triplet [U, L, C] = [1, 0, dmn], 
m =2, , M. 

4. Create a return arc (la , lb) with capacity cost triplet [U, L, C] = [1, 1, 0] to 
transform the network G = (V, E). 

A circulation network representation for an M-machine clustering problem is shown 
in Fig. 4.1. The flow on arc (roa, mb) is always 1 so that each node mb can supply 1 
unit of flow to one destination node na (n * m), by the flow conservation condition. 
The network solution is equivalent to minimiTing the total sum of dissimilarity values 
between the machines used in the clustering problem There is a formation of one main 
loop which is a closed loop from the source node (lb) to the sink node (la). 
Additionally, sub-loop(s) is (are) also formed vriiich correq)onds to a closed chain not 
containing both the source node and the sink node. These closed loops correspond to 
the machine cells. 

4.1.2.1b Circulation network construction for generalized grouping 

The distance parameter or the negative similarity parameter defined in section 4.1.1 is 
used as the per-unit cost of flow and the corresponding network flow problem then 
becomes the minimum cost flow problem (MCFP). 

In order to construct the circulation network for process family formation in 
generalized grouping problem following procedure is proposed. 

1. Create two nodes represented by Pka and Pkb for each part k,k = 1, ,K. 

Also, arbitrarily let node Plb be the source node and Pla the sink node 

of the resulting network G = (V, E). 

2. Create two nodes represented by u and ib for each process plan i = 1, , q. 



Rest of the network is structured according to the following rules. Each arc is 
assigned with capacity-cost triplet values [U, L, C] to indicate an upper bound 
on its flow, a lower bound on its flow and the per unit flow cost, respectively. 

3. For each node Pka, k 1, create an arc directed from Pka to Pkb with capacity- 
cost triplet defined by [U, L, C] = [1, 1, 0]. 

4. From each node Pkb, create arcs directed to each ib, iiFk with a capacity-cost 
triplet [U, L, C] = [1, 0, 0] . 

5. For each nodeia,iIFk, k?^!, create an arc directed fi'omiato ib with 
capacity cost triplet [U, L, C] = [1, 0, 0] . 

6. For each of the nodes ib, iiFk , k = 1, ,K create (q - |Fk|) arcs directed 

firom ib to each node ja, jlFk, with a capacity-cost triplet [U, L, C] = [1, 0, dy]. 

7. For each of the nodes ib, create an arc directed fi'om ib to ja with capacity- 

cost triplet [U, L, C] = [0, 0, B], iJIFk, k = 1, , K where ‘B’ is a very 

large positive number; or, equivalently, do not create any arc directed firom 
ibtoja, iJeFk. 

8. For each node ia, ieFi, create an arc directed to the sink node Pka with 
capacity-cost triplet [U, L, C] = [1, 0, 0]. 

9. Create a return arc (Pka, Pkb) with capacity-cost triplet [U, L, C] = [1, 1, 0] to 
transform the network G = (V, E) into a circulation network. 

A circulation network representation for a K-part and q-process plans clustering 
problem is shown in Fig. 4.2. The arcs (Pkb, Pka) has the capacity-cost triplet [U, L, C] 
= [1, 1, 0] which indicate that the flow through these arcs are always 1 so that each 
node Pkb can supply 1 unit of flow on one of the arcs (Pkb, ia), ieFk, therby the node ib 
connected to ia by arc (ia, ib) can supply 1 miit of flow to the destination node ja, 
jgFk, by the flow conservation condition. The cost of shipping 1 unit of flow firom 
each node ib, ieFk to each destination node ja, j^Fkis the dissimilarity coefficient dy. 
Thus the network solution is equivalent to roinimiTing the total sum of distances 
between the plans involved in clustering along with selecting only one plan of each 
part. 
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4.1.2.1c Equivalence of the proposed model (MAAM) to the UCCNF 
model 

Let the flow ja) on arc (L, ja), ieFk, j^Fk be represented by Xy and the 
flow f(ia, ib) on arc (L, L) by x;;. Then the objective function of the proposed model 
which minimizes the total dissimilarity between the plans transforms into minimum cost 
flow problem for the circulation network constructed in the previous section. The 
transformed objective function is 
minimize 


z 

ieF, 

k 


Z dij f (ib Ja) 


z 


l^j IJ D 


z 

leF, 


d..f(i 

n 


a' V 


Since h(^ib, ja) = 0 and di; = 0, finally the objective function is 
minimize 

The constraints are as follows: 

(I) constraint (4.4) becomes 


(4.8) 


Z Z = 1; for all part k= 1, ,K 


(n) constraint (4.5) becomes 


(4.9) 


Z Z ia) = 1 


(4.10) 


(HI) constraint (4.6) becomes 


Z Z f(ib ^ ia) ^ ^ ^ ^^Jb ' ia - 0 or 2 


i eFk, j ^Fk for all nodes representing process plan. 


(4.11) 



Studying the network reveals that each node h, ieFk can either supply 2 imits of flow 
or 0 mut of flow depending upon whether its corresponding ib is supplied 1 imit of flow 
or not from the node Pkb. If it gets the supply, then corresponding ia will supply 1 imit 
of flow on arc (ia, ib) so that ib in turn can supply 1 unit of flow (constraint 4.9) to its 
destination node ja, j^Fk and the other unit of flow to node Pka for which it receives 
unit supply from node jb (constraint 4.10) and becomes jb’s destination node. Thus the 
constraints (4.4) , (4.5) and (4.6) together indicate the flow conservation conditions 
for the nodes of the network. Constraint (4.7) is for lower and upper bound on the 
flow. 

Solution of the network results into formation of closed loops. Among these one will 
be the main loop containing both the nodes ib and ia of plan ieFi which are connected 
with source node Plb and the sink node Pla respectively as (Plb, ib) and (ia, Pla). Rest 
of the loops termed as sub-loops do not contain both the source and the sink node. 
These loops corresponds to the process families. 

4.1.2.2 Solution methodology 

The model formulated through equation (4.1) to (4.5) is a 0-1 integer LP problem 
which can be solved by any procedure solving 0-1 integer LP model e.g. LESTDO or 
CPLEX. Solution of the model results into some non-zero variables which is then 
interpreted as fo rmin g some loops. The number of loops gives the number of process 
family formed. 

We give here an illustration. The table below shows the parts and their process plans. 
Each part ma y have more than one process plan. The resulting solution is in terms of 
the reference number of the plan. 

Example problem 

The problem with 5 parts (K = 5), a total of 1 1 process plans (q = 1 1) and 4 machines 
(M = 4) is diown ia the Table 4.1. 



Table 4. 1 : Input Data for the Exainple Problem 


part 

k 


1 

1 

1 

2 

2 

3 

3 

4 

4 

5 
5 


plan 
number 
of part 
k 


1 

2 

3 

1 

2 

1 

2 

1 

2 

1 

2 


plan machines 

reference 
number 

1 

12 3 

1 0 0 1 

2 0 10 

3 11:0 

4 011 

5 10 1 

6 1 0 0 ' 

7 0 10 

8 10 0 

9 1 ' 0 1_ 

10 0 0 1 

11 ,10 0 , 


4 

1 

1 

0 

0 

0 

1 

1 

1 

0 

1 

0 


The distance coefficient matrix D = [dy] calculated for the above problem and used 
the objective function is shown in Table 4.2. 



Table 4.2 : Distance matrix for the example problem 


plans 1 2 3 

1 0 00 00 

2 00 0 00 

3 00 00 0 

4 ^ 1 2 2 

5 2 4 2 

6 2 2 2 

7 2 0 2 

8 2 2 2 

9 2 2 4 

10 0 2 4 

11 3 3 1 


4 5 6 7 

2 2 2 2 

2 4 2 0 

2 2 2 2 

0 00 4 2 i 

00 0 2 4 

4 2 0 ' 00 

2 4 - 00 0 

4 2 2 0 

2 0 2 4 

2 2 2 2 

3 113 


8 9 10 11 

2 _ 2 ’ 0 3 

2 4 2 3 

2 2 _ 4 1 

4 2,2 _ 3 

2 0 2 1 

0 2 2 1 

2 4 2 3 

0 00 2 1 

00 0 2 1 

2 2 0 00 

1 1 00 0 


After solving the problem using the proposed formulations, MAAM (eqs. 4.3 - 4.7) by 
CPLEX we get the following resuhs. 

X27 = 1, X72 = 1, X59 = 1, X911 = 1, X1105 = 1; and all other Xy’s are zero. 

From this we can infer that plan reference number 2 and 7 form a loop and plan 
nmnber 5, 9, and 11 form another loop (see Fig 4.3). Thus two process families (2,7) 
and {5,9,11} are formed. The result exactly matches with that by Kusiak’s model 
Correspondingly part families are; {1,3} and {2,4,5}. 

4.2 IMPROVEMENTS IN THE KUSIAK’S FORMULATIONS 

Referring to the Kusiak’s p-median model [1987] as discussed in section 3.1 the 
constraints of assigning process plans to a process family (eqn 3.7) is 


xij < xjj 


foraIli,j 


( 4 . 12 ) 






Through this constraint Kusiak postulated that at least one plan should be selected for 
the process family formation. This implies that even one process plan can make a 
process family. But it is not a very practical assumption from the point of view of GT. 
concept. In fact, in GT problem we do not usually consider and include that process 
plan which is disjoint of all other plans since there is no similarity of the manufacturing 
attributes of this plan with any other plans. Hence, it is argued in agreement with the 
GT concept that at least two process plans of two dififerent parts should be selected 
out of the several plans available for a process family formation. Accordingly, the 
above constraint is modified as; 

Ex =0; ifxjj = 0. 

1=1 

>1; ifxj„=l. (4.13) 


The integer programming formulation for the above constraint in terms of inequalities 
is given as. 



q 

E X.. < Bx.. . 

i=i y jj 

for all i, j 

(4.14) 

and 

E X.. > 2x.. . 
i=l iJ JJ 

for all i, j 

(4.15) 


where B is the ma ximum number of parts desired to be kept in a family. 

The constraints (3.6) of the Kusiak’ s model for the formation of ‘p’ required number 
of process families is 

q 

Z xjj < p. (4.16) 

j=l 

The advantages of the constraints (4. 14) and (4. 15) in place of constraints (3.8) are of 
three fold. 

1. This modification will not require the constraints (4.16) for the formation of a 
prespecified required number of process families. Thus it results into natural 


grouping. 



2. If one wishes to use the constraint (4.16) when the number of process femilies is 
given in advance it can be included as well Thus it provides the flexibility between 
imposed grouping and natural grouping. 

3. Besides these the number of the constraints is also decreased in comparison to that 
of Kusiak’s model. 

4.2.1 Improved Models 

Based on the above improvements and observations, the following two models are 
proposed. 

(1) Imposed Grouping model (IMPGRP) 

This results by replacing the constraints (3.7) by the constraints (4.14) and (4.15) in 
the Kusiak’s model (section 3.1). The complete formulation is presented as 
The objective of the model is to form the process plan families by maximiTing the 
value of similarity measure among them. 

q q 

Maximize sij xij (4-17) 

subject to 

(I) Selection of only one process plan for a part 

Z Z xii = 1 (4.18) 

ieFkj=l 

(n) Formation of maximum required number of groups 

Z xjj < p (4.19) 

j=l 


(IQ) AssignmCTit of process plans to groups 



q 

Z X.. < Bx- 

i=l u JJ 

for alii, j 

(4.20) 

and 

q 

Z X.. > 2x.. 

i=l u 

for all i, j 

(4.21) 


(IV) Binary decision variables 

jQj = 0 or 1 (4-22) 



( 2 ) Natural Grouping model (NATGRP) 

This results by eliminating the constraints (3.6) and replacing the constraints (3.8) by 
the constraints (4.14) and (4.15) in the Kusiak’s model (section 3 . 1 ). The complete 
fomulation is presented as : 

The objective of the model is to form the process plan families by maximi^Tng the 
value of similarity measure among them. 

q q 

Maximi 7 .e Z Zgiixii (4-23) 

i=lj=l •* •’ 

subject to 

(I) Selection of only one process plan for a part 

.Z Zxii = 1 (4.24) 

(H) Assignment of process plans to groups 

q 



Z X.. < Bx.. 
i=l iJ JJ 

q 

for all i, j 

(4.25) 

and 

Z X.. > 2x.. 

i=i y a 

for all i, j 

(4.26) 


(m) Binary decision variables 


Xij = 0 or 1 (4.27) 

A comparison between the models will follow later in the chapter. 

4.2.1. 1 Solution methodology 

The iiBproved models EMPGRP and NATGRP are still an integer programming 
formulation \\iiich can be solved by any procedure solving 0-1 LP model e.g. LINDO 
or CPLEX. Solution of the problem results into some non-zero variables which are 
interpreted for femily formation. 

Example problem 

Take the exairple of Table 4. 1 and use the distance matrix D. 

The result we get in both the cases of natural grouping (NATGRP) and inqrosed 
grouping (IMPGRP) with required process femihes, p = 2 , are the same and are given 
below. 

X 27 = 1 , X 77 = 1 , X 59 = 1 , X 99 = 1 , Xiio 9 = 1 and all other jqj’s = 0 . 



The process families are then {2,7} and {7,9,11} and part femilies are {1,3} and 
{2,4,5}. These results exactly match with the Kusiak’s solution. 

Some example problems ( section 7.3) are presented in chapter 7 where the advantages 
of improved models EMPGRP and NATGRP can be clearly seen. 

4.3 TYPES OF GENERALIZED GROUPING PROBLEMS 
Problem Type I : Part-Wise Routeing Flexibility (PWRF) 

This type of problem is generally of small order size ( <10 parts).In this type, more 
than one type of operation of a part can be carried out either on a single machine or on 
more than one machine. This means that the set of alternative machines for a particular 
type of operation of the manufacturing system is different for different parts. This type 
of problem is generally of small order size (<10 parts). Exanq)le problem given in 
Table 4. 1 is an example of this type of problem. 

Problem type II : Operation-Wise Routeing Flexibility (OWRF) 

This type of problem is for bigger manufecturing environment. In this type, alternative 
machines are available system-wise, that is the set of alternative machines for a 
particular type of operation is the same for all the parts over the whole manufacturing 
system. Example problem given in Table 7.7 is an example of this type of problem. 

4.4 A H EURISTIC PROCEDURE 

The proposed model, MAAM, (section 4.1) and the natural grouping model, 
NATGRP, (section 4.2. 1) may result into a large number of process femilies for larger 
problem. So a heuristic method is applied for mer^g some of the process femilies 
vsdiich uses the same machines or the difference is only due to the alternative machines. 
Later the assignment of machines to the family is done while targeting the objective of 
minimizmg intercell movements. A refinement procedure is also included which 
reassigns bottleneck parts to further reduce the intercell movement. The procedure 
consists of the following steps. 



step 1 Identification of sets of ahemative machines. 

This step is useful only for generalized grouping problem of type II (OWRF). 

When multiple routeing of part is possible because of presence of alternative machines 
for a type of operation. Identify those machines for different type of operations and 
make their sets. . 

step 2 Merging the part families. 

If the machines used by a process family is a subset or superset of the machines used 
by any other process family, merge the two process families. If subset or superset 
criteria fails between the two process Emilies only due to the presence of machines of 
the same set of alternative machines then merge them by assigning either of the two 
machines. Restructure the process family accordingly. Repeat for every other pair of 
process family. 

step 3 Assigning the machines to the process fennlies. 

(i) Compute for each machnie m the usage factor Umi for each process family 1. 


(ii) Assign the machine to the process family ‘1’ wfiere Umi is maximum, the ties are 
broken arbitrarily. 

special case : Generalized grouping problem of type n (OWRF) will involve the 
following step also. 

(iii) If Sumi for an alternative machine is zero i.e. no operation has been assigned to 
it, then place it in place of that machine of the same set which has Umi > 0 for more than 
one process family It may be noted that such a situation arises vAen number of 
alternative machmes for a type of operation is relatively large. 

step 4 Merging the process families and machine cells. 

Merge two process families and their correq)onding machine cells if it eliminates 
intercell movement and does not violate the cell size limit imposed. Stop when no 
merging is possible or the groups formed are digoint. 
step 5 Refinement procedime 


Go to step 3. Assign the bottleneck parts to that cell which has maximum number of its 


operations. Ties are broken arbitrarily. 


weNT 




k A. J 


Note : 


Simple grouping and generalized grouping type I problems need not to go through 
stepl and step (iii) step 3. 

Kusiak and Agarwal’s model need to go through steps 3 and 5. 

4.5 COMPARISON OF MATHEMATICAL MODELS 

A comparison of Kusiak’s model (PMM), Agarwal’s modified model, the proposed 
model (MAAM), and the proposed in:q)rovements in kusiak’s model (IMPGRP and 
NATGRP) can be made on the basis of niimber of constraints, number of variables and 
the order of complexity. Order of con^jlexity, as such, is calculated for the basic 
operations of an algorithm which are used repetitively. In that sense we do not have 
basic operations in an integer programming formulation as optimization is a 
simultaneous process. However, the number of aU possible combinations which are 
likely to be explored while solving the problem through different integer programming 
formulations can be a good yardstick of measuring the conoplexity of the formulation. 
This has been discussed in the following paragraph. 

The number of ways in vdiich n different objects can be partitioned into p non-empty 

subsets is given by Stirling’s number, S(n,p) = p” / p! (Duda and Hart, 1973). In most 

grouping problems, exploration of the total number of all feasible grouping alternatives 

for the optimal solution is sinply not possible, even when large coiiq)uters are 

available. The number of feasible clustering alternatives may be astronomically large 

even for relatively small grouping problems. For exan:5)le, the total number of possible 

ways of grouping 20 objects in five groups is approximately 794,728,600,000. The 

combmatorial nature of the problem becomes more important when the number of 

groups is unknown, vdiich is so m most groupmg problems. In fact, this is what is to be 

determined m most of the problems. It can be seen that m this case the total number of 

n 

possibilities is given by HS{n,j). For the 20-object exan: 5 )le, it would be of the 

7 = 1 



order of 10*^. A corr^dete enirmeration for the optimal solution should therefore be 
ruled out. This highly combinatorial nature of the grouping problem restricts the use of 
Kusiak’s model to relatively small size problems. The proposed model applies the 
concept of associating one object with the other without having the number of clusters 
known. In that case the total enumeration of all possible combinations for n objects 
n(n-l). For 20-object example , this number would be 180. Thus, it can be concluded 
that the proposed model (MAAM) is highly con 5 >utationaIly efficient even for the large 
scale problem Similar discussion can be made for the natural grouping model 
(NATGRP) also. 

Following tables shows comparison between the models explicitly: 

Table 4.3 : Total nmaber of enumerations for various models 

I - 

Models 

Kusiak’s model 

(PMM) 

Agarwal’s modified model 

Proposed model 

(MAAM) 

Imposed grouping model 

(IMPGRP) 

Natural grouping model 

(NATGRP) 
where 

K 

n = total number of ways selecting one process plan for each part k. 

k = \ 

and K Cp Total number of combination of selecting p elements out of K 


Total number of enumerations 
K 

n I Fk !• S (K, p) 
k=l 

K 

n iFkl*Kcp»S(K-p,p) 
k = l 

n |Fkl*K(K-l) 
k = l 

K 

n |Fkl*Kcp*S(K-p,p) 
k = l 

n iFki»K^ 

k = l 


elements (K > P). 



Table 4.4 compares the models on other criteria. 


Table 4.4 : Comparison of various models 


Criteria 

Kusiak’s 

model 

(PMM) 

Agarwal’s 

modified 

model 

IMPGRP 

model 

NATGRP 

model 

Proposed 

model 

(MAAM) 

GT concept 

Violates 

Violates 

Violates 

Agrees 

Agrees 

Kind of 
^ouping 

Imposed 

grouping 

Inq)osed 

grouping 

Imposed 

grouping 

Natural 
: grouping 

Natural 

groupmg 

Number of 
variables 


; 

q2 


q2 

Number of 
constraints 

q^+K+l 

q+K+1 

2q+K+l 

2q+K 

2K+q 

Heuristic 

application 

steps 1,3 
and 5 

steps 1,3 
and 5 

steps 1,3 
and 5 

All steps 

All steps 



5 

SIMn.ARITY AND DISTANCE MEASURES 


In this chapter we discuss some similarity and distance measure developed along 
with some existing one from literature. These measures has been developed in the 
need of various types of inputs like the traditional binary input matrix, process time 
input matrix, operation sequence of process plans along with machine utUizaticm 
penalty inputs. 

We have observed for the various techniques that it is the measure of similarity (Sy) or 
dissimilarity (dy) between the process plans vriiich plays crucial role in optimization. So 
better or stronger is this measure, more practical the result will be. Also for different 
types of inputs e.g. binary input, process-time input, routeing sequence input and 
machine utilization input etc., one should have some measure which will be able to 
rightly account for the similarity or the dissunilarity between the process plans. We 
propose three types of coefficients for different types of inputs described in the 
following sections along with their parallel coefficients available in the literature. 

5.1 MACH INE WEIGHTED SIMILARITY COEFFICIENT 
(MWSC) 

This is defined for the binary matrix input. A part either requires a machine or does 
not, Le. 

aim ~ indicates plan i uses machine m; 

= 0, otherwise. 

The weightage to a ma chin e is given by the sum of I's under that machine in the 
incidence matrix. The significance of this is that the machine \\frich is required by 
larger number of plans is more likely to be a general purpose machine on which a big 
investment has occurred. So, obviously the management would like to utilize it highly. 



(5.1) 


The weightage for a machine m Wm is thus given by 

N 

^m ~ Sami 
i = l 

for any two plans i and j 

N 

MWSCjj — Z!wm-^(aim»ajm) 
i = l 

= -oo; IJ eFk 

= 0;ifi=j (5.2) 

where 5(auii^ aj^) = 1 ; ^ aim = ajm ^ ^ 

= 0 ; otherwise 

5.2 PROCESS TIME RELATED COEFFICIENTS 

Till now we considered 0-1 part machine matrix input to calculate the coefBcient 
between two plans. Considering process time of operation as input, one should have 
some defined coefiBcient. ha the following paragraphs we propose some coefiBcients 
based on the process time. 

5.2.1 Multiplicative Weighted Similarity Coefficient (MuWSC) 

Mosier (1985) defined three types of similarity coefiBcient between machines which use 
process time of operation as input. In parallel to his mult^licative wei^ted similarity 
coefiBcient between the pair of machines we propose a similarity coefiBcient between 
the pair of process plans. 

The proposed coefiBcient is defined for operation process times as inputs. The 
similarity between a pair of process plans is the difference of the weights of aU the 
common machines required by them and weights of the machines required by any one 
of the two process plans. The measure evaluates the similarity between the plans by 
taking products of the weight of the common machines and subtracting the weights of 
the machmes required by any of the two plans. 

Let tim = total processing time of plan i on machine m (tin, > 0 only wBien aim = 1). 



For the confutation of the proposed similarity coefiScient MuWSCij, weight of 
machine is taken as the relative weightage on the basis of its requirements with respect 
to overall requirements of machines and is given as 

q 

2 tim 

w = i=l 

“ NNON. / 

where NNONim = the total number of non zero elements in the matrix. 

The numerator of Wjj^ represents the total processing time on machine m over all 
plans. The denominator denotes the average time load over all machines. 

The Multiplicative Weighed Similarity Coefficient between a pair of process plans is 
given by 


MuWSCij 


MM2 

2 tim-t jm ~ £ (tii)i~tj]ii) 

m = 1 m = l 

M 

swii 

m=l 


(5.4) 


The &st term hi the numerator of the MuWSCjj is a measure of the requirement of the 
same machines by plans i and j. similarly the second term of the numerator is a measure 
of the misses of the plans with regard to machines. The term in the denominator is to 
standardize the coefficient in the range of [0, 1]. MuWSCy can be used in place of Sjj in 
all the objective functions for maximization in the formulations of chapter 4. 

5.2.2 Time of Movement Distance Coefficient (TMDC) 

This is the sum of times of those operations of plan i and plan j vriiich are not 
common. 

Let the function tjm) = tj^; if tim = 0 

= 0; otherwise 

This function gives the time of intercell movement when the machines required for the 
plan j is not available in plan i 



Let TOMmax ~ the maxi mum total time of movement \\4ien operations of plan j are 


carried out on the machines corresponding to plan L 

M 

= Max[ S f(aim»aiiii) ] 

L j m = 1 


(5.5) 


The distance coefficient TMDCij, when the operations of plan j are carried out on the 
machines corresponding to plan i, is given by 


M 


2 f(aiiii,a jm) 

TMDC;; = “ = ^ 

TOM 


(5.6) 


max 

It can now easily be inferred from the above definition that 

TMDCij * ™DCji 

It is because machines requirement by the two plans i and j may differ. 

Apphcation of TMDC will result into assignment of machines to a process family to 
minirnize the time of interceU traffic. TMDCjj can also be used in place of dy in aU the 
objective fimction for the minimization in the formulations of chapter 4. 

5.2.3 Weighted Part and Time Similarity Coefficient (WPTSC) 


Application of MuWSC does ensure the minimization of intercell movement of parts as 
it tries to keep those parts together which have similarity of types of operations in 
terms of machine requirement. But this may not ^ultaneously ensure that the time of 
stay outside the cell due to this intercell movements are also a minimum. Hence, the 
TMDC should also be taken mto consideration to have a tradeoff between the 
minirnization of number of part movements and the minimization of time of movement. 
This is to be noted that we are only considering the time of operation while 
determining the time of movement not the time of transportation between the cells. 
Hence, another coefficient, to be called as Part and Time Similarity Coefficient 

(WPTSC), can be defined as follows. 

WPTSCij = Ip * MuWSCij - It * TMDCij (5-7) 

where Ip and L are the m 5 )ortance factors for part and time of movement, respectively. 
Clearly, WPTSQj WPTSCji 



5.3 OPERATION SEQUENCE BASED COEFFICIENTS 

Operation sequence is considered to be one of the important criteria to judge the 
similarity between the process plans. This has an inqiortant influence on the cell layout 
and thereby on the material handling cost and time. The process plans having the 
highest similarity between them when grouped together will enhance the flowability of 
the components and will reduce the load on the problem of scheduling. Before 
presenting the proposed coefficient we discuss some of the earlier proposed 
coefficients grabbing the similarity between the process plans due to operation 
sequence. 

5.3.1 Transformation Based Similarity Coefficients (TBSC) 

(Tam, 1990] 

Tam (1990) described three types of transformations viz. substitution, insertion, and 
deletion to capture the operation sequence similarity between the process routeing of 
parts. The operation sequence of a plan is represented by a string having the indices of 
machines required by a part. The similarity coefficient defined by Tam specifies the 
similarity between the parts as the similarity between these strings m terms of the 
minimum number of weighted transformation required to carryout the conversion. 

Definitions 

Let the conversion of plan i to plan j require at least the following transformations: 

ns : number of interchanges( substitutions) in the string of plan j 

n; : number ofmsertions to the string of part j 

n<i : number of deletion firom the string of part j 

Let Ml, M 2 and M 3 be the wei^ts assigned to ns, ni and xu respectively. 

The similarity coefficient between plan i and plan] belonging to dififerent parts is: 

Sij = l/min(Mins + M 2 ni + Many) (5.8) 



5.3.2 Sequence Length Based Similarity Coefficient (SLBSC) 
[Choobineh, 1988] 


Here it is assumed that similarity between part routeing is proportional to the length of 
common operations sequences between them. Hi^er the number of common 
operations required by conqionents in sequence, higher is the similarity between them 
Definitions 

The SLBSC is defined as the sum of similarity between the common operation 
sequences of length 1 to Le between any two parts. 

Let 

Le : the maximum length of sequence of operation considered for comparison 
Xoi = 1, if operation o is performed on plan i; 

= 0, otherwise. 

Cij (y) ; number of common operation sequences of length y between plan i and plan j’s 
routeings. 

Li : length of routeing of plan i 

Lnin : mini (L) = number of operations in the smallest part routeing 

The similarity is defined as 

M 


Sij 


J_ 

Le 



Z Xoi Xoi 
0=1 


Z (xoi + Xoi“Xoi Xoi) 
0 = 1 


Le 

+ z 

y = 2L 


Cjj(y) j 

min-y + 1 


(5.9) 


The first term in the similarity coefihcient defined is an another way of representing 
Jaccard similarity coeflScient which is defined as the ratio of number of machines 
required by both plans to the sum of the machines required by both and the machines 
required by either of them So this term represents the similarity of common operations 
between the plans. In efiTect, it represents the similarity due to single leagth of 
operation sequence. The second term calculates the similarity between the common 
operation sequences of the length 2 to Le. This similarity coeflScient is dependent upon 
the choice of L*. Its high value while results in more accmrate measurement of 
similarity coefihcient, it involves more con^utational time. Generally its value is kept 



equal to the maxi r o u in length of common sequence that exist between plans of two 
different parts. 

5.3.3 Operation Sequence and Machine Utilization based 
Dissimilarity Coefficient (OSMUDC) 

The proposed OSMUDC is as follows. A process plan has certain sequence of 

operation which is called process routing sequence. It is not necessaiy that two process 

plans which visit the same type of machines may have the same sequence of visiting 

them. The similarity coefficients discussed above do not consider the sequence 

possibility in the reverse direction vdiich is likely, however, to avoid this we can put a 

high penalty factor on reverse flow of a part in the cell layout. In the case of SLBSC, 

the Jaccard coefficient does take the reverse flow into account in some sense but not 

the similarity of sequence in reverse direction. Sometime it is found that if a few back 

flow in the middle of flow of the component through the layout is allowed, then it 

might be possible that afterwards the flow is smooth in forward direction which is not 

accounted by TBSC and SLBSC. Hence, a more practical similarity from the 

manufacturing point of view can he assessed by the similarity of the sequence of 

operation in both forward as weU as backward direction, however, to again avoid 

backward flow we penalize the back flow of part. To account the similarity due to 

similar types of ma chin es used by two routeings, we allow reverse flow of coirqionent. 

The explanation of terms used for defining the dissimilarity coefficient is given below. 

Definitions 

Intercellular moves 

Given the layout of machines of a machine cell, a part moves within the cell according 
to its requirement fi)r operation in a certain sequence. These moves from machine to 
machine within the cell are termed as intracellular moves. 



Intracellular moves 

Given the layout of machines of a machine cell, a part may have certain requirement of 
machines for which the machine is not available in the cell and for that it has to visit 
some other machine cell. Every such move is termed as intercellular move. 

Utilization penalty junction 

This is assigned to a machine by the management and depends on how and to what 
extent the management wants the machine to be utilized under dijBferent demand 
scenario for its products. A general purpose machine, on v^ch higher investment is 
incurred diould be maximum utilized. Similarly a fecial pmpose machine should be 
aimed at to be idle for the least period. When the demand stracture of products 
changes, a readjustment is supposed to be performed. The processing capability and 
flexibility of a machine should always be kept in mind while determining this. It is quite 
possible that the process plans also get redefined in such cases. In this sense we can say 
that this involves complete design of GT problem It becomes a complex function 
when factors hke cost, life, maintenance etc. are also accoimted besides the above 
mentioned complexity of the environment. In general we can assume the utilization 
factor to be equal to the fi'equency of occurrences of the machine in the ^stem 
Whenever, two or more machines are deared to be kept together in one cell their 
utilization factor should be equally high. 

5.3.3.1 Calculation of moves, penalty and dissimilarity coefficient 

Notation 

Wb = penalty fector for back flow ; Wb ^ 1. 

Wj = wei^tage for one intracellular move. 

Wq = weightage for one intercellular move. 

Ei = set of machines required by plan i 
UPjjj = utilization penalty for machine (m = 1 to M) 


TUPij = ITJPm, meEi and m^Ej. 



Table 5.1: Intracell and intercell movements for plan j as per plan i 


Machine for 
planj 

Intracell movements 

latercell 

raovement 

Remarks 

1 

4 

no 

forward flow 

2 

3 

no 

back flow 

3 

1 

no 

forward flow 

4 

1 

no 

- do - 

5 

- 

yes 


8 

3 

no 

forward flow 

7 

- 

yes 



Thus the total intraceU movements, Xij = 4 + 3.Wb +1 + 1+3 

Another way of counting intracell movement can he to count ‘0’ movement for the 

forward movements and actual movements for the backward movements. That is, 

Xij = 0 + 3 + 0 + 0 + 0 = 3. 

It can be noted that in either case Xjj Xjj 

While going through layout we also observe that locations of machines 5 and machine 
7 of plan j are not foimd. So for those machines we count as number of intercell 
movement. That is the number of machines of plan j which is absent in plan i for the 

present case is Yy = 2 . 

Also, Yy ^ Yjj 

This also can be observed that machines 6 and 9 of plan i are not being utilized by the 
plan j so we charge utilization penalty for them Here, the sum of the penalty of those 
machines of plan i which are not utilized by plan j, TUPjj = UPg + UP 9 . It can be 

easily seen that TUPjj ^ TUPjj. 



The dissimilarity coeflBcieot, d'ij, \\iieii plan j proceeds in the machine layout of plan i 


wi-Xij+woYij ^ 
wi-Xjuax Wo-Yjjjgx 


TUP.. 

y 

Z UP, 
ma. 


m 


(5.10) 


Again, it can be noted that d'ij ^ d'ji 

Hence the net dissimilarity between plan i and plan j is, 

dij = dji = (d'ij + d'ji)/2 (5.11) 

From the final result of grouping we wiU have a sound indication of making the layout 
of machine for all process family. For making the final layout of machine cells we 
follow the following rule. 

For a process family determine for its every plan i 


Di = Z d'ij (5.12) 

j i 

where i and j belong to the given family. 

The layout of ma chin es for the process family will be according to that plan wiiose Dj 


is minimum of all. 
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MEASURES OF EVALUATION 


Practical problems of group technology are usually and unfortunately framed according 
to the requirement of the methodologies. Many of the aspects of production (e.g. part 
machine incidence, their process time, operation sequence, multiple routeings, machine 
capacity constraint, production rate, demand, availability of material handling 
equipment, and many other such related to a manufacturing environment) are 
discarded. However, in some cases inq)ortant ones are tried to get accounted for while 
developing a techmque. Consideration of all those aspects which have an impact on the 
final solution will end up with a too con 5 )lex problem to even formulate. So, at least 
some of them which play a greater role in decision making should be considered for the 
purpose of framing the problem Hence, solution of practical problem without 
adequate generalization and theoretic formulation is supposed to be useful only in 
specific situation for which the problem has been framed. Such results can not usually 
be apphed to other situation with different problem scenario. Group technology has 
suffered from this malady for a long time. When researchers and practitioners claim 
success with a certain method it is necessary to evaluate the result on an absolute 
quantitative scale. The practical use of such a measure is that subjective claim can be 
compared objectively. 

The measures selected here evaluate the solutions of binary matrix, process time input 
matrix, and part routeing sequence. The measures to evaluate the binary matrix tends 
to compare the solution with perfect block diagonalized matrix and those for routeiag 
sequence coixpare with a hypothetical solution with zero intercell traffic. For 
evaluating the solution of process time input, and the operation sequ^ce input 
global evaluation criteria as suggested by Harhalakis et aL(1990 a) will be used. 



Another measure for solution of process time input can be the number of parts 
movement along with the time of interaction among the cells. 


6.1 EFFICIENCY MEASURES TO EVALUATE THE 
SOLUTION CONSIDERING MACHINE PART 
MACHINE MATRIX 

The measures to compare the part machine matrix take a real value according to the 
row and column arrangement of part machine matrix. The output matrix provided by 
cell formation algorithm happens to be arranged in a way that ah the non zero elements 
are concentrated towards the diagonal of the matrix. These efficiency measures 
evaluate the concentration of non zero elements (I’s) in the matrix. The more Ts are 
arranged along the diagonals the better is the output depicted by the value of efficiency 
measures. The various measures incorporated in this work are grouping efficiency 
[Chandrasekharan and Rajagopalan, 1985], grouping efficacy [Kumar and 
Chandrasekharan, 1990] and relative efficiency [Chandrasekharan and Rajagopalan, 
1987]. These are reproduced synoptically in the following subsections for the purpose 
of evaluation of various grouping approaches discussed in the Chapter 4. 

6.1.1 Grouping Efficiency 

The concept of grouping efficiency is developed to provide a quantitative standard on 
a rational scale for comparing diffe rent solutions by different approaches to the same 
problem. 

Let A = [aij] be the part-machine binary (0-1) incidence matrix, i = 1, ,N and j = 

1, M. 

Let {Dj}, vshere r = 1,2 R, be the diagonal submatrices obtained by associating 

the product group with the machine group in the diagonalized matrix (R — total 
number of such groups). A perfect grouping iirq)Ues that 
aij = 1 if ajj e (Dj) and 

= 0 if aij ^ {Dr) r = 1 R 



The diametrically opposite situation, therefore irq)lies that 
ajj = 0 if ajj e {Dj.} and 

= 1 ifajj ^ {Dj-} r= 1 R 


The eflSciency scale is defined to cover these extreme cases, its value being set at one 
for the former case and zero for the later. 

Goodness of grouping depends on two aspects: within group utilization and intercell 
movement. From the matrix point of view the concentration of non zero elements in 
the diagonal submatrix refers to utilization and the presence of such elements outside 
the diagonal submatrices represents intercell movements. A better grouping therefore 
means an increase in the utilization, a decrease in the intercell movement or both. The 
number of non zero elements in the diagonal blocks is given by 


Mr 

s 


R 

e^ 2 2 ^ 2 ^ ajy 

r = li = Mr-l + lj = Nr-i + l 


Nr 

s 


(Mo = 0,No= 0 ,Mr = N,Nr = M) 

Nj- = Sum of number of parts cluster 1 through cluster r. 

Mj- = Sum of number of machines cluster 1 through cluster r. 
and the number of non zero elements outside the diagonal block is 


(6.1) 


“"O 


= e - ea 


Where 'e' is the total number of non-zero elements (I’s) in the matrix. 

and the number of voids in the diagonal blocks is 

R 

ev = ZMrNr ' 
r = l 

Grouping efficiency can therefore be considered as a weighted average of two 
efficiencies Tj j and t| 2 defined as follows. 


ni = 


ed 

R 

( SMfNr) 

r = l 


( 6 . 2 ) 



n2= 1- 


(6.3) 


eo 

R 

(MN- ZMfNr) 
r = l 


The numerator of the expression for t] 1 is the number of non zero elements in the 
diagonal block and the denominator the total number of elements there in. Similarly ri2 
is the ratio of the number of zeros in the off diagonal block to the total number of 
elements therein. 

The grouping efficiency 't|' can be expressed as the wei^ted average of ti j and 7)2 

Ti = qrii + (l-q)Ti2 (6.4) 

T| satisfies the basic requirements of non dimensionally, non - negativity and the range 
of 0 to 1. The weighing factor enables the analyst to alter the enphasis between 
utilization and mtercell moves. A high value of q ensures greater utilization. 
Limitations of grouping efficiency 
1 . Effects of weights 

The weighted average of rn and x\2 gives the analyst aiough freedom to decide the 
relative importance between intercell movements and voids in diagonal blocks. 
Chandrasekharan and Rajagopalan [1985] concluded that a value of q = 0.5 leads to 
the situation of equal weights to voids and exceptional elements. An analysis of the 
expression would reveal that this is not tme. When large matrices are diagonalized, for 
any number of blocks greater than two it can be observed that the first term in 
expression has a much smaller denominator conq)ared to the second term, whereas the 
numerator are more or less of the same order. As the matrix size increase this disparity 
between the first and the second term widens and for large and / or sparse matrices the 
second term becomes less and less effective. In fact if q and 1-q are of the same order, 
the effect of intercell movement is never reflected in the efficiency value in the case of 



large and sparse matrices. Therefore for equal weight to voids and exceptional 
elements it is necessary to choose low value q. Further it is desirable to bring the two 
parts of the expression under a common denominator. 

2. Effect of sparseness of matrix 

Choice of q would be more rational if it is linked with the size and qjarsity of the 

matrix. As a result, q is chosen as the density of an imag inary matrix that has an 

identical and perfect block diagonal structure as the given matrix. 

K 

q = EMrNr/MN (6.5) 

r = l 

The efficiency value reduces to 

T) = 1 - (e^ + Cq) / MN (6.6) 

This simplified expression assumes that the perfect diagonal form has 100 percent 
efficiency which is reduced by the presence of voids and exceptional elements. 

The foEowing characteristics are inqjortant ui the case of the modified expression for Tj 
after choosing q as discussed earlier. 

The main difference between this ejqpression and the earlier expression for efficiency is 
that the present expression has the size of the matrix (MN) as its denominator, 
whereas in earlier expression tj 1 and 112 have two different denominators. It is also 
interesting to note that it is possible to assign relative wei^ts for voids and 
exceptional elements, because it is in the present form that voids and exceptional 
elements have equal wei^ts. 

A major defect of grouping efficiency in both fiirms is its low discriminating power. It 
is obvious that (e^ + Cq) is much less than the size of the matrix (MN). Thus for large 
matrices the term (ey + e^) / MN wiE take values of the order of 0.00 1 which leads to 
an efficiency close to one. 

Another serious defect Ees in the definition of zero efficiency. Zero has been defined as 
the most degenerate case where (ey + Cq) = MN. This nrpEes that the diagonal 
blocks are fuB of zeroes and off diagonal and off diagonal are fuE of ones. It can be 
seen that this definition is realistic firomphyacal point of view. 



3. Effect of the size of the matrix 

By accepting e as the denominator of the eJBBciency function it will have the 
form 

r) — 1 - [(Oy + eQ)/e] 

= (®d - ©vVe 

The main drawbacks of this function are as follows : 

1. Zero point of the eflSciency value is not properly defined. The value of "n becomes 

0 when Cjj = Cy i.e. 

R 

ed = ey = l/2[( ZMrNr)] 
f = 1 

This means that zero efliciency is dictated by the feet that there are 50 percent 
voids in the diagonal blocks irrespective of the number of exceptional elements 
present. 

2. A major defect of groupuig efficiency (both in earher and modified form) is its 
discriminatmg power. It is obvious that (ed + ey) is much less than the size of the 
matrix. Thus for a large matrices the term (ey + eo)/MN will take values of the 
order of 0.001 which leads to an efficiency close to one. 

3. This function takes negative values when ey<ed. 

6.1.2 Grouping Efficacy 

For a large denommator the discriminating power of the efficiency function is poor and 
for a small denominator the function may itself become negative. 

The meaningful denominator between 'e' and 'MN* is the operational zone of the 
matrix. The operational zone of block diagonal form consists of the diagonal blocks 
usually consists of both ones and voids, and the exceptional elements consists of 
only ones. 

Operational zone = (Voids + ones) in the diagonal blocks + ones in off diagonal 
blocks. 

= ey + ed + Cq 

Oy "t" e 



Thus expression can be modified as 

”n ~ 1 ■ [(Cy Cq) / (Cy + e)] 

— (e - e^) / (ey + e) 

This is the grouping efficacy (GEF). 

GEF = (1 -(p)/ (l+4>) (6.7) 

(p = Number of exceptional elements / total number of operations 
(j) = Niunberofvoids in diagonal blocks /total number of operations 

6.1.3 Relative Efficiency 

The limiting efficiency states that for any block diagonal stmcture defined by sequence 
of row clusters and column clusters, there exists a limit on the maximum attainable 
efficiency. 

Let Mj and Nj- be the number (dimension) of vectors in the rth row cluster and cohinro 
cluster respectively. 

Arrange the sequence of integers Tf and 'M* so that 

Mi^M2^ Mr 

and 

Ni^N2 ^ Nr 

Then accommodation 'S' in the diagonal block is maximum when M and N are 

correlated to form a diagonal submatrix. Then 

R 

^max ~ SMj-Nr 
r = l 

The accommodation 'S' is minimum when sequence of 'M and 'N' are reversed and the 

corresponding cluster M and N are correlated to form diagonal blocks. Then 

R 

S min ” SMi-Nr 

r = l 

for any general combination 

Smin - S ^ Sma x 

The grouping efficiency (Chandrasekharan and Rajagopalan, 1985) is ^ven by 

Ti = qrii + (l-q)T12 



Casel ; (e<Sinin) 

In this case voids will appear in blocks, but exceptional elements can be hypothetically 
avoided. Hence the hmiting value of eflSciency suffers because of first component of 
grouping efficiency y\i. 


Case2 : (OSmax) 

In this case exceptional elements are inevitable, althou^ voids can be theoretically 

avoided. Hence the second component of grouping efficiency t \2 suffers. 

Tio = q + [(l-q)*(MN - e/MN - 
Cases : ( S - < S < Sjjjax ) 


mm 


For any value of ‘S’ between S min and SxnaY h is theoretically possible to find a 
sequence so that 


S»e 


In this case, Tio = 1. 

Relative efficiency is defined as the ratio of grouping efficiency and limiting efficiency. 
This ratio is a useful check to find out whether a given algorithm has achieved the best 
of what is achievable under the constraints of given data. 


6.2 GLOBAL EFFICIENCY MEASURES 

This class of measures evaluate the result vdiich are either in the form of binary matrix, 
process tim e matrix or in the form of part routing sequence. Three measures wduch 
have been designed to evaluate their results, by Harhalakis et. aL(1990a) are as 
follows. 

6.2.1 Global Efficiency 

This is the ratio of the total number of part operation that are performed within their 

respective cells to the total number of operations in the system. 

N 

Zsi 

Global efficiency = 



where sj is the number of operations in the rj that are performed in the cell 
corresponding to part i. 

Global efficiency examines the grouping power of any method in a absolute way. It 
reflects the effectiveness of the assignment in confining the operations of parts within 
their respective cells. It assigns a penalty for each operation to be performed on a part 
outside the cell. 

6.2.2 Group Efficiency 

This criterion evaluate the efficiency of the solution to minimize the intercell moves. 
The intercell moves are maximum, when every machine belongs to a single group and 
are zero when aH machines are in a single group. 

The ratio of ihe difference between the total nmnber of maximum external cells that 
could be visited and the total number of external cells actually visited by all parts to the 
total number of maximum cells that could be visited by them is known as group 
efficiency. 

The maximum number of external cells that could be visited is 

N 

Ew= ZMin{(qi-l),(w-l)} (6.9) 

i = l 

where qj is the number of different machines belonging to the same routing rj 
Xjc = 1 if part i visits cell c 

= 0 otherwise c = 1 w 

The total number of cells actually visited by all parts, 

N w 

A„= Z Z(Xic-l) (6.10) 

i= lc = 1 

Group efficiency = (‘ 'D 

Jbw 

This ratio assigns the same penalty whether there is one or more operations that are 
performed on the part in the same foreign cell This is a more abstract measure, but it 



reflects the effectiveness of the assigntnent to confine the parts to as few foreign cells 
as possible. 

In this measure of evaluation the penalty given account for the number of cells in 
which the external operations are performed. Global efficiency gives the penalty 
regardless of the number of external cells visited. So the higher the group efficiency, 
the better the assignment is. There must be at least two cells in the system for the 
group efficiency to be calculated. 

6.6 Group Technology Efficiency 

It is the ratio of the difference between the ma ximuTn number of intercell traffics 

possible and the number of intercell traffics actually being required by the system to the 

maximum number of intercell traffics possible. 

The maximum number of intercell traffics possible in the system is 
N 

1= ECri-l) (6.12) 

i = l 

Xjo = 0, if operations o and o+l are performed in the same cell; 

= 1, otherwise. 

where o = 1 rj - 1, for part i. 

The number of intercell traffic actually required by the system is 
N ri-l 

u= z 2 Xio 
i=lo=l 

Hence, Group technology eflSciency = (I - U ) / 1 (6. 14) 

This is a very powerful criterion, which takes into account the sequence in which the 
operations are performed apart from the cell in which they are performed. This 
efficiency provides a penalty of 1 for each intercell move. Alternatively, this efficiency 
can be evaluated by using the actual production volume of the parts (uj) wiule 
computing the intercell moves, I and U, thus providing a weighed penalty, it can be 
viewed as a measure reflecting the saving in intercell moves. 
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IMPLEMENTATION, EXPERIMENTATION, 
RESULTS, AND OBSERVATIONS 

This chapter presents its contents in three phases. 

Phase I : E?qpeiimentation and comparison of the proposed model (MAAM) with 
Kusiak’s and Agarwal’s modified models. 

Phase n : Validation of the proposed similarity and dissimilarity coefficients. 

Phase in : Testing of the proposed improvements in Kusiak’s model (PMM). 

To cany out three-phase experimentation, an integrated system is developed and 

implemented on HP 9000 machines in C language. The system gives mathematical 

formulation for any combination of the model and similarity/dissimilaxity coefficient 

desired. All the five models (Kusiak’s PMM, Agrawal’s modified model, MAAM, 

IMPGRP and NATGRP) have been included. Four coefficients viz. Lee’s dissimilarity 

coefficient, machine weighted similarity coefficient (MWSC), weighted part and time 

similarity coefficient (WPTSC), and the operation sequence machine utilization 

dissimilarity coefficient (OSMUDC) are incorporated in the system. Thus, in total 5x4 

= 20 different types of possible fiirmulations are experimented with. The ^stem is 

menu driven in which user can select the combination of the model and the coefficient 

to be used to formulate the problem. The formulation thus resulting is stored in an 

output file which is solved by CPLEX. Computational CPU times are recorded. The 

output fi-om CPLEX is interpreted to form the process families. The proposed heuristic 



is then applied manually to merge the process families and then to foim part families 
and machine groups. Various grouping efficiency measures for the final solution are 
computed manually. 

• Phase I makes use of the integrated system for formulating some example problems 
using the proposed model (MAAM), Kusiak’s model (PMM) and Agarwal’s 
modified model. These models are then compared on the criteria of CPU-time 
performance and efficiency measures. 

• Phase n validates the proposed sirnilarity/dissiinilarity coefficients by testing some 
example problems containing the inputs required for the calculation of the 
coefficients. The testing makes use of the system to apply either the proposed 
model (MAAM) or Kusiak’s model (PMM) for sohing the problem. 

• Phase in makes use of the system for testing some exanple problems on IMPGRP 
and NATGRP models. The results are than compared with those obtained by 
applying original Kusiak’s model (PMM). The basis of comparison is again CPU- 
time and efficiency measures discussed in Chapter 6. Some noteworthy 
observations are also discussed. 

PHASE I 

7.1 COMPARISON OF PROPOSED MODEL (MAAM) WITH 
KUSIAK’S AND AGARWAL’S MODIFIED MODELS 

The exanple problems taken are binary (0-1) plan and/or part machine incidence 

matrix. The problems are divided into two categories: (a) ample grouping problem 

with each part having single process plan, and (b) generalized groping problem with 

parts having alternate process plans. 



7.1.1 Simple Grouping Problem 

The 0- 1 part machine incidence matrix for jSve sunple grouping exanqjle problems are 


shown in Tables 7. 1 to 7.5. 

Example Problem 1 

Table 7. 1 : Part-Machine Incidence Matrix for Example Problem 1 

Machines 

123456789. 10 

1 001 101 0000 

2 1000001001 

3 0 100100 100 

40001010010 
5 0 100100 100 

P6 001001 0010 

a 7 0 00000 1 00 1 

rSOlOOlOOlOO 
t 9 0 0 1 1 0 1 0 0 1 0 

s 10 100000 1001 

11 1 0 0 0 0 0 1 0 0 1 

12 100000 1001 

13 0 1 0 0 1 0 0 1 0 0 

14 001 101 0010 

15 0 1 00 1 0 0 1 00 


Example Problem 2 

Table 7.2 : Part-Machine Incidence Matrix for Exan^le Problem 2 

Machines 

1234567891 11111 

0 1 2 3 4 5 

I'O 10001 10001 100 0 

2 000010001100101 
P 3 101001000000010 

a 4 101101000000010 

r 5 0 00000011100101 

t 6 010001100011000 

s 7 0100010000110 0 0 

8 1011010 0 0010010 

9 001101100000010 

10 0 000000 1 1 1 0 01 0 1 



Example Problem 3 

Table 7.3 : Part-Machine Incidence Matrix for Example Problem 3 

Machines 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

1 











0 

1 

1 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

2 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

3 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

4 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

1 

5 

0 

1 

1 

0 

0 

1 

0 

1 

0 

1 

0 

6 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

7 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

8 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

9 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

1 

10 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

11 

1 

1 

0 

1 

0 

0 

0 

0 

0 

1 

0 

12 

0 

1 

1 

0 

0 

1 

0 

1 

0 

1 

0 

13 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

15 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

16 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

17 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

1 

18 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

19 

0 

1 

1 

0 

0 

1 

0 

1 

0 

0 

0 

20 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

21 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

22 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 



ExaiTiple Problem 4 

Table 7.4 ; Part-Machine Incidence Matrix for Exan:q)le Phroblem 4 


1 

1 0 
2 0 

3 0 

4 0 

5 0 

6 1 

7 1 

8 0 

P 9 0 

a 10 0 

r 11 0 

t 12 0 

s 13 0 

14 0 

15 0 

16 0 

17 0 

18 0 

19 0 

20 0 
21 0 
22 0 

23 0 

24 0 


Machines 
2 3 4 5 6 7 8 

0 0 110 10 
0 0 1 0 0 1 0 

1 1 0 0 0 0 0 

1 1 0 0 0 0 0 

0 0 0 0 0 0 1 

0 0 0 0 0 1 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 1 0 1 

0 0 0 0 1 0 1 

0 0 0 0 1 0 0 

0 0 0 0 1 0 1 

0 0 0 0 0 0 0 

0 0 0 0 1 0 1 

0 0 0 0 1 0 1 

0 0 0 0 1 0 1 

0 0 110 10 
0 0 0 0 0 0 0 

0 0 1 0 0 1 0 

0 0 110 10 
0 1 0 0 0 0 0 

0 0 0 0 1 0 1 

0 0 1 1 0 0 0 

0 0 0 0 0 0 0 


9 11111 
0 12 3 4 

0 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 1 0 0 0 

1 0 0 0 0 0 

0 0 0 0 1 0 

0 0 0 1 1 0 

0 0 0 1 1 0 

1 0 0 0 0 1 

1 0 0 0 0 0 

0 0 0 0 0 1 

1 0 0 0 0 0 

1 0 0 0 0 1 

0 0 0 0 0 0 

1 0 0 0 0 1 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 1 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 1 0 0 0 

0 0 0 0 0 0 

0 0 0 0 1 0 

0 110 0 0 



Example Problem 5 

Table 7.5 : Part-Machine Incidence Matrix for Example Problem 5 


Machines 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

1 

1 

1 

1 

1 

1 











0 

1 

2 

3 

4 

5 

6 

1 

0 

0 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

2 

0 

1 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

3 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

5 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

6 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

7 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

8 

0 

0 

0 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

10 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

11 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

17 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

18 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

19 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

20 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

21 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

22 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

23 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

24 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

1 

1 

0 

0 

0 

25 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

1 

0 

0 

0 

0 

28 

0 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

32 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

33 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

34 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

35 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

36 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

37 

1 

1 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

38 

0 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

39 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

40 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

41 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

42 

1 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

43 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 



7.1.1.1 Methodology, solution and observations 

We find the solution of above problems using Kusiak’s PMM model, Agarwal’s 
modified model and the proposed model (MAAM). The coefiBcients used are: (a) Lee’s 
distance coefficient, and (b) machine weighted similarity coe£Bcient(MWSC), 

For the proposed model number of part famihes to be formed is not prespecified. The 
proposed heuristic procedure is appUed (excluding stepl, \yhich is appUcable with 
generalized grouping). 

For Kusiak’s and Agarwal’s modified models the number of ma xiTmim number of part 
famihes (or cells) is assumed and solution gives the part families formed. Hence, to 
reach at the best solution possible we run the model for a few times (trial and error). 
Step 3 of the heuristic proposed is used for machine assignments to part families. The 
refinement procedure of the heuristic improves the solution \\iienever it is possible to 
improve. 

Solutio_ns of Example Problem 1 

Problem size : (15x10), total number of operations (e) = 46. 

Solutions are shown in the subsequent tables (Table 7.6 and 7.7). 



(a) Lee’s distance coefficient 


Table 7.6 ; Solutions of various models using Lee’s distance coefficient for Example problem 1 


Model 

P 

CPU 

Part 

Machine 

Number of 

Number 

Remarks 


number 

time 

f^ilies 

cells 

bottleneck 

of 



of loops 

(sec) 



operations 

voids 



1 




(eo) 

(ev) 


Kusiak’s 

p = 4 

0.06 

(1,4.6, 

(3,4,6,9) 

2 

3 

Family of a single 

model 



9,14) 




plan is formed. 

(PMM) 



(2,10, 

(1,7,10) ■ 







11,12) 








(7) 

(-) 







(3,5, 

(2,5,8) 





- 

-- 

8,13,15) 






p = 3 

0.05 

(1,4,6, 

(3,4,6,9) 

0 

4 

No single plan 




9,14) 




family. 




(2,7,10, 

(1,7,10) 







11, 12) 








(3,5,8, 

(2,5,8) 




- 



13, 15) 





Agarwal’s 

p = 4 

3.06 

(1) 

(-) 

3 

3 

Single plan 

modified 



(4,6,9,14) 

(3,4,6,9) 



family is formed. 

model 



(2,7,10, ■ 

(1,7,10) 







11 12) 








(3,5,8, 

(2,5,8) 





- 


13,15) 






p = 3 

1.55 

(1,4,6, 

(3,4,6,9) 

0 

4 

No single plan 




9,14) 




family. 




(2,7,10, 

(1,7,10) 







11, 12) 








(3,5,8, 

(2,5,8) 







13, 15) 





Proposed 

1 = 6 

0.01 

- do- 

- do- 

0 

4 

Cell size limit = 5 


model 


(MAAM) 

Observations 

1. There is remarkable difference of coiBputational time of three formulations. The 
proposed model (MAAM) takes the least time (0.01 sec.) and Agarwal’s model 
takes the ma viTTmm time 3.06 sec. with four families (p ~ 4) and 1.55 sec. with 3 


families (p = 3). 



2. can observe the GT concept violation by both Kusiak’s and Agarwal’s modified 
models with p — 4 as single plan famihes (7) and (1) respectively are formed. 

3. The sum total of the number of bottleneck operations (Cq) and the number of voids 
(Ov), that is, (Oo + Ov) is 6 for Agarwal’s model (p = 4), 5 for Kusiak’s model (p = 
4), and 4 for the remaining cases of the Table 7.1. Based on (e^ + e^) the solutions 
obtained by Kusiak’s and Agarwal’s modified models both with p = 4 is worst as 
(go + ev) is high. 

4. Bottleneck operations are more critical than void elements from the point of view 
of practical result. Hence, based on the number of critical operations (eo), the 
solution obtained by Kusiak’s and Agarwal’s modified models, both with p = 4 are 
worst as eo = 2 and 3 respectively for them, whereas for the proposed model 
(MAAM) there is no bottleneck operation. 

5. Thus, in the case of Kusiak’s and Agarwal’s modified models with higher value of 
(say p = 4), while the values of (eo + Cv) are higher i.e. the grouping quahty is bad, 
simultaneously the famUies of single process plan also get formed and 
computational time also increases. 

6. Only merging step of the proposed heuristic is reqtiired to get the final solution of 
the proposed model (MAAM). 

7. The solution by the proposed model exactly matches with the best solution 
obtained by Kusiak’s and Agarwal’s modified models both for p = 3. 

8. It can be argued that the total CPU time taken by the proposed model and heuristic 
together will still be lower because aheady the problem is reduced to loops which 
is at least half of the number of parts and the merging step have the order of 

complexity as 0(1^). 



Efficiency measures for the best solutions obtained by Kusiak’s and Agarwal’s 

modified models both with p = 3, and the proposed MAAM model are; 

grouping efficiency 97.33% 

grouping efficacy 80.76% 

global efficiency 100.00% 

group efficiency 100.00% 


(b) MWSC 

By using MWSC in the three models we achieve the best result (eo + ev) as obtained 
for p = 3 in Kusiak’s and Agarwal’s modified models and the proposed model in (a). 
Even using higher value of p (say p = 4) in Kusiak’s and Agarwal’s modified models, 
solution results into forming three part families. The CPU time is, however, different as 
shown in Table 7.7. 

Table 7.7 : Computational times of various models using MWSC for example problem 1 


Model 

p or 1 

CPU time ( 

Kusiak’s model 

p=4 

0.03 


p = 3 

0.05 

Agarwal’s modified model 

p = 4 

36.22 


p = 3 

20.30 

ftoposed model 

1 = 4 

0.01 


Observations 

1. In comparison to (a) lesser CPU time for Kusiak’s model and higher time for 
Agarwal’s modified model are needed to solve the problem 

2. For Kusiak’s and Agarwal’s models both, no family of single process plan is 
formed with even hi^er value of p (say p = 4 or 5) and solution results into such a 
number of process fa mili es (less than p) wbich gives the best result possible (Le. 
with p =3). 



3. Because of this, CPU time for Kusiak’s model with higher value of p (say p = 4) is 
less than that required for the best result possible i.e. with p = 3. 

4. Less number of loops are formed for the proposed modeL It is because of the 
somewhat uniform distribution of operations on machines in the example problem. 
That is most of the machine’s weighing foctor lies in a small range. For the large 
range case, i.e. large variability of operations assignments on machines, larger 
number of loops are expected. 

5. Why to use MWSC ? 

• Generally it is found that interceU moves are more due to the machine having more 
number of operations, so the aim is to keep those parts together creating these 
moves due to these machines. 

6. Where to not use MWSC? 

• When a problem has its total operations on machines distributed in a large range. 
That is one machine has few operations while some other machine has large 
number of operation assigned. 

Solutions of Example Problem 2 

Problem size : (10x15) and total nmnber of operations (e) — 48. 

Solution are shown in the Tables 7.8 and 7.9. 



(a) Lee’s distance coefficient 


Table 7.8 ; Solutions of various models using Lee’s distance for example problem 2 


Mcxiel 

porl 

CPU 

time 

(sec) 

Part 

families 

Machine 

cells 

Co 

©v 

Remarks 

Kusiak’s 

model 

(PMM) 

p = 3 

0.02 

(1,6,7) 

(2,6,7, 

11, 12) 

1 

10 

No single 
family plan 

(2,5,10) 

(5,8,9,10, 

13,15) 

(3,5,8, 

13, 15) 

(1,3, 4,6, 
7,14) 



Aganval’s 

modified 

model 

p = 3 

0.2 

- do- 

-do- 

-do- 

- do- 

-do- 



Proposed 

model 

(MAAM) 

1 = 3 

0.02 

- do- 

-do- 

- do- 

- do- 

-do- 


Efficiency measures are: 


grouping efficiency 

92.66% 

grouping efficacy 

81.03% 

global efficiency 

97.91% 

group efficiency 

95.00% 


(b) MWSC 

We get the same results obtained in Table 7.8 but the computational times are different 
given in Table 7.9. 

Table 7.9 : Computational time of various models using MWSC for example problem 2 


Model 

p orl 

CPU time (sec) 

Kusiak’s model (PMM) 

p = 4 

0.03 


P = 3 

0.02 

AgarwaTs modified model 

p = 4 

1.98 

P = 3 

0.86 

Proposed Model (MAAM) 

1 =3 

0.01 


Solutions of Example Problem 3 

Problem size : (22x11) with totalnumber of operations (e) = 78. 
Solutions are shown in the following tables (Table 7. 10 and 7. 1 1) 





(a) Lee’s distance coefficient 

Table 7. 10 ; Solutions of various models using Lee’s distance coefficient for example prdjlem 3 


Model 

p or 1 

CPU 

Part 

Machine 

Go 


Remarks 



time 

families 

cells 






(sec) 






Kusiak’s 

p = 3 

0.29 

(1A3.7,11, 

(1.4,5,10) 

15 

10 

Refinement procedure 

model 



15,16, 




could not improve Cv 

(PMM) 



20,21,22) 




with p = 3 


or 

0.16 

‘ (46,9,10, 

(7,8,9,11) 



because the procedure 




14,17,18) 




does not consider the 




(5,8,12, 

(2,3,6) 



void minimization 


p=5 


13, 19) 




which can, however, be 
added. 

for p = 5 merging step 
and refinement step 
reduces Cv and eoboth. 

AgarwaFs 

p = 3 

< = 30 

- do- 

- do- 

“ do^ 

- do - 

- dc - 

modified 

or 







model 

p = 5 

A 

!l 

'o 

- do- 

-do- 

- do- 

- do- 


Proposed 

1 = 9 

6.03 

(1,2,3,7,11, 

-do- 

14 

10 

Cv is reduced by 

model 



15,16,20, 




applying refinement 

(MAAM) 



21,22) 




procedure. 




' (4,9,10,14, 
17,18) 
(5,6,8,12, 




Cell size limit = 5 




13,19) 






(b) MWSC 


Table 7.11: Solutions of various models using MWSC for example problem 3 


Model 

p or 1 

CPU 

Part 

Machine 

eo 


Remarks 



time 

families 

cells 






(sec) 






Kusiak’s 

p = 3 

' 0.15 

(1,2,3,7,11, 

(1A5,10) 

15 

10 

Refinement 

model 



15,16,20, 




procedure could not 

(PMM) 

or 


21,22) 




improve result as 



0.13 

(46,9,10, 

(7,8,9,11) 



it does not consider 


p = 5 


14,17,18) 




void minimization 




(5,8,12, 

(2,3,6,) 



which can be 




13, 19) 




added. 

Agarwal’s 

P = 3 

A 

II 

U) 

o 

- do- 

-do- 

-do- 

-do- 

- do- 

modified 

or 







model 

p = 4 







Proposed 

1 = 8 

0.03 

(1,2,3,7,11, 

- do - 

14 

10 

ev is reduced by 

model 



15,16,20, 




applying the 

(MAAM) 



21,22) 




refinement 




■ (4,9,10,14, 




procedure. 




17,18) 








■ (5,6,8,12, 




Cell size limit = 5 




13,19) 







Efficiency measures for (a) and (b) both are : 


Efficiency measures Kusiak’s and Proposed model 


modified (MAAM) 

Kusiak’s model 

grouping efficiency 89.67% 90.08% 

grouping efficacy 73.11% 73.91% 

global efficiency 87.18% 87.18% 

group efficiency 73.68% 73.68% 


Observations 

In aU the above three problem it is observed that 

1. In most of the cases, CPU times with MWSC are found to be less than then- 
corresponding cases of Lee’s distance. In case of the proposed model less number 
of loops are formed for MWSC and simultaneously CPU time is also less which is 
due to somewhat uniform distribution of operations on various machines. 

2. The proposed model takes minimum computational time as well as the final 
' solution (Go = 14, Co + ev = 24) is better than that by Kusiak’s and Agarwal’s 

modified models (eo = 15, eo ev = 25). 

3. In Agarwal’s modified model we had either large computational time or had to 
limit this. 

Solutions of Example Problem 4 

Problem size : (24x14), total number of operations (e) = 60. 

Solutions are shown in Tables 7.12 and 7.13. 



(a) Lee’s distance coefficient 


Table 7. 12 . Solutions of various models using Lee’s distance for example problem 4 


Model 

por 1 

CPU 

Part families 

Machine 

Co ' 

ev 

Remarks 



time 


cells 






(sec) 






Kusiak’s 

p = 4 

0.38 

(1,2,17,19, 

(4,5,7,10) 

29 

3 

' Refinement does not 

model 



20,23) 




improve Cv and Co 

(PMM) 



(3,4,6,7,8,18, 

(2,3,11) 



further. 




21,24) 








(6,7,8,18) 

(1,12,13) 







(5,9,10,11,12,13, 

(6,8,9,14) ' 




- 



14,15,16,22) 


- 



Agarwal’s 

p = 4 

< = 

- do- 

-do- 

29 

3 

with refinement 

modified 


50 





procedure otherwise 

model 












- 

- 


\ e^ = 30, Co = 4. 

Proposed 

1 = 9 

0.07 

- do- 

- do - 

- do- 

-do 

Refinement 

model 






- 

prcK^edure reduces Cv 

(MAAM) 







further. 


- 


- 




Cell size limit =5 

Efficiency 

measures are; 







EflSciency measures All models 


EflSciency measures All models 
grouping eflSciency 90.47% 

grouping eflScacy 64.44% 

global eflSciency 95.08% 

group eflSciency 91.89% 


It was observed that Kusiak’s model (p =3) results into three families (with eo — 53 and 
ev = 3) which is worse than the result with p = 4 in Kusiak’s model (with eo = 29 and e\ 
= 3, see Table 7.12). 


(b) MWSC 

Table 7.13: Solutions of various models using MWSC for example problem 4 


Model 

porl 

CPU 

time 

Part 

femilies 

Machine Co 

cells 

ty Remarks 

Kusiak’s 

model 

(PMM) 

p = 4 

(sec) 

0.23 

(1,2,17,19, 

20,23) 

(3,4,6,7,8, 

18,21,24) 

(5,9,10,11, 

12,13,14, 

15,16,22) 

(4,5,7,10) 53 

(1,2,3,11, 

12,13) 

(6,8,9,14) 

3 After refinement 

procedure 



Aganval’s P ~ 4 < = 50 same as (a) 

modified 

model 

Proposed 1=11 0.05 - do - 

model 

(MAAM) 


same as(a) 29 3 No improvement 

by refinement 
procedure 

- do - - do - - do - Refinement 

procedure reduces 
ev and Co both. 
Cell size limit = 5 


Observations 

In this problem also, use of MWSC with these models results in less CPU time. 
However larger number of loops results for the proposed MAAM model. Number of 
loops resulting from the proposed model is a direct consequence of the way the 
distribution of operations on machines. Larger number of loops results when the range 
of distribution of operations is large or the distribution is skewed, that is, the difference 
between minimum number of operation and the maximum number of operations 
assigned to the machines in the system is large. The measure of large variability of. 
operations distribution on various machines can be determined by the difference 
between the number of operations assigned on various machines and the average 
nxunber of operations on a machine in the system, which is the ratio of the total number 
of operations (e) and the total number of machines (M). 

EjB&ciency measures are 


Efficiency measures Proposed and Kusiak’s model 

Modified model 

grouping eJBficiency 90.47% 83.33% 

grouping efficacy 64.44% 50.42% 

global efficiency 95.08% 95.08% 

group efficiency 91.89% 91.89% 


In one go. Only the proposed model (MAAM) and the proposed heuristic together 
gives the op timal number of part families for wduch the efficiency measures are the best 



possible. Kusiak’s model has to be used few times to reach at this result. This can be 
validated by the reduction of eflBciency measures with p = 3 and using Kusiak’s model 
(see above table). 

Solutions of Example Problem 5 

Problem size : (43x16), total number of operations (e) = 126. 

Solutions are shown in the subsequent tables (Table 7.14 and 7.15). 

(a) Lee’s distance coefficient 

Table 7. 14 : Solutions of various models using Lee’s distance for example problem 5 


CF*U Part families 

time 

(sec) 

3.30 (1,6,7,12,13,17,25,26, 

31,34,35,36,39) 

' (2,4,10,18,28,32, 

37.38.40.42) 
(3,11,20,24,27,30) 

(5,15,16,22, 

29.33.41.43) 
(8,9,14,19,21,23) 

3.60 (1,6,7,12,13,17,25,26, 

31,34,35,36,39) 
(2,4,10,18,28, 
32,37,38,40,42) 
(3,11,20,24,27,30) 
(5,8,9,14;i5,16,19,21, 

22,23,29,33,41,43) 


Agarwal’s p = 7 

< = 

- do- 

modified 

50 


model 

<= 

same as Kusiak’s 

p = 5 

50 

p=5 


Proposed 1=18 0.18 

(1,12,13,25,26,31,39) 

model 

(2,4,10,18,28,32, 

(MAAM) 

37,38,40,42) 


(3,11,20,24,27,30) 


(5,8,9,14,15,16,19,21, 

22,23,29,33,41,43) 

(6,7,17,34,35,36) 


Machine cells 

eo 

e.. 

Remarks 

(3,6,7,10,14) 

39 

74 

After 

refinement 

(1,2,9,6) 



procedure 

otherwise 

(8,11,12,13) 



even worse 

(5,15) 



result. 

(4) 




(3,6,7,10, 14) 

29 

79 

After 

refinement. 

(1,2,9,6) 



procedure. 

Without 

(11,12,13) 



refinement 

(8,4,5,15) 



step Co = 38 
andev=84. 

-do- 

- do - 

- do 

merging step 



- 

and refinement 




procedure. 

sameas 

39 

74 

After 

Kusiak’s 



refinement 

1! 

L/t 



procedure. 

(7,10) 

30 

58 

Best result is 

(1,2,9,6) 



achievable. 

Cell size limit 

(11,12,13) ■ 



. =7. 

(4,5,6,8,15) 



; Big dimension 
problem is 

(3,4) 



solved more, 
efficiently 


Model p or 1 


Kusiak’s p = 5 

model 

(PMM) 


p = 4 



Efficiency measures are 


Efficiency measures 


grouping efficiency 
grouping efficacy 
global efficiency 
group efficiency 


Kusiak’s and 
Modified model 
Cv = 79, Co = 29 
84’30% 
47.31% 
76.98% 
61.33% 


Proposed model 
(MAAM) 

87.21% 

52.17% 

76.19% 

57.14% 


(b) MWSC 

T able 7.15: Solutions of various models using MWSC for example problem 5 


Model 

p or 1 

CPU 

Part families 

Machine cells 

eo 

Cv 

Remarks 



time 








(sec) 






Kusiak’s 

p = 5 

16.92 

(1,6,7,12,13, ' 

(3,6,7,10,14) 

29 

79 

After refinement 

model 



17,25,26,31, 




procedure 

(PMM) 



34,35,36,39) 

(2,4,10,18, 

(1,2,9,6) 



Without refinement 




28,32,37, 




step Co = 38 andev = 




38,40,42) 

(3,11,20, 

24,27,30) 

(11,12,13) 



4 

1 




(5,8,9,14,15, 

16,19,21,22, 

(8,4,5,15) 







23,29,33, 








41,43) 





Agarwal’s 

p = 5 

<= 

no integer 





modified 


50 

solution 





model 



possible. 





Proposed 

1 = 20 

0.17 

(1,12,13,25, 

(7,10) 

30 

58 

Best result is 

model 



26,31,39) 




achievable. 

(MAAM) 



(2,4,10,18, 

(1,2,9,6) 



Cell size limit = 7 




28,32,37, 

38,40,42) 

(3,11,20, 

(11.12,13) 



Big dimension 




24,27,30) 




problem is solved 




(5,8,9,14,15, 

16,19,21,22, 

(4,5,6,8,15)' 



more efi5ciently. 



- 

23,29,33, 

41, 43) 
(6,7,17,34, 

(3,4) 







35,36) 






Efficiency measures are the same as that for case (a). 



♦ In case (a), Kusiak’s model with p = 4 gives a better result than with p =5. In case 
(b), with p = 5 we get the same result of Kusiak’s model with p = 4 of case (a). In 
the case (b) with p —5 in Kusiak s model the resulting solution forms four families. 

♦ Proposed model forms five families and gives the best result as the value of eo + 
Cv (= 88) is smaller than that for Kusiak’s model with p = 4 (that is its better 
solution) which is 108. 

♦ CPU time is very less as compared to that of Kusiak’s model 

7.1. 1.2 General observations 

1. Where do we need to go through refinement procedure? 

Mostly with MWSC we need this procedure, however, if the distribution of 
operations is not sparse tlien we may not need this. It is because of the fact that 
high weigliing factor of the machine having larger number of operations when 
maximized will pull that part towards itself ( that is towards the cell containing this 
machine) which uses some other machines of low weighing fectors besides the high 
weighing factor machines. So, even when the part has more number of operations 
in some other cell it will be found, as a result, in the cell having the machine of hi^ 
weighing factor. 

2. It is because of the refinement procedure appHed to reassign some leftover 
bottleneck parts that final result obtained by the three models are quite near to each 
other. 

3. Use of MWSC results in less than the number of families ‘p’ given wften p is given 
a higher value whereas with Lee’s distance coefficient the solutions results into 


exactly ‘p’ number of families given. 



4. Use of MWSC with the proposed model results in less number of loops when 
distribution of operations on various machines is somewhat uniform. 

5. Use of MWSC in case of sunple grouping can be advantageous when the 
distribution of operations on machines is less sparse. Otherwise, one needs to go 
through the refinement procedure arises. 

6. Bearing the above observations in mind, use of MWSC in the case of generalized 
grouping is not advisable since when there is more than two alternative machine for 
the same operation it may result in having Sumi > 0 for only two of these 
alternative machines, where 1 is the family to which machine m has been assigned. 
Also Zumi for one of these two wfll be higher than the other thus creating 
unnecessary intercell movement which, otherwise, by using Lee’s distance 


coeflScient would not occur. 



7.1.1.3 Graphical analysis of solutions 

Graph of number of plans and ‘e’ vs CPU time for various models is shown m Fig. 7.1. 


FIGURE 7.1 iNumber of parts( or e) vs CPU time 

for simple grouping 

4 - ; — □ — Kusiak I 



10 15 22 24 43 

48 46 78 60 126 

Number of parts ( or e) 


It is found that for the proposed model (MAAM) CPU time varies quadratically with 
number of plans while for Kusiak’ s model (PMM) it varies e5q)onentially. No such 
conclusion can be drawn for Agarwal’s modified model as we had to limit its time of 
computation for most of the example problems taken. 

Graphs between problem number and CPU time for Kusiak s model and the proposed 
model using Lee’s distance coefficient and MWSC are shown in Figs. 7.2 and 7.3. 
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FIGURE 7.3 : Problem number vs CPU time 
(Proposed MAAM model) 
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It is observed that for Kusiak’s model time variation is large vshen MWSC results in 
less number of process families than ‘p’, the maximum number of famihes required. 
However, for the proposed model this variation is neghgible. 

7.1.2 Generalized Grouping Problems 

Five example problems have been taken in the class of generalized grouping problems 
(2 of PWRF type and 3 of OWRF type). The 0-1 input matrix for the first three 
exanqjle problems (2 of PWRF type and 1 of OWRF type) are shown in the following 
tables (Tables 4. 1, 7. 16 and 7. 17). 

Example Problem 6 ( problem type PWRF ) 

The input data are the same as in Table 4.1. 

Exauiple Pipblem 7 .( problem tiq)e PWRF ) 

Table 7. 16 : Binary Input Incidence Matrix for Exaoqjle Problem 7 


3 4 5 


1 1 0 
1 0 0 
1 0 0 
1 1 0 
1 0 0 
1 1 1 
0 1 1 
1 1 1 
0 0 1 
1 0 0 
1 0 0 
1 0 0 
1 0 0 
1 0 0 
0 1 0 


Part 

k 

1 

1 

1 

2 

2 

3 

3 

4 

4 ' 

5 

5 

6 
6 
7 
7 


Plan Plan 

no. of refer, 
jrait k no. i 


Machines 
1 2 


1 

2 

3 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


1 

1 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 


1 

1 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 



Example Problem 8 ( problem type OWRF ) 


Table 7. 17 : Binary Input Incidence Matrix for Exanq)le Problem 8 

Machines 


Part 

plan 

Plan 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

k 

no. 

refer. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 


of 

no. 






















part 

i 






















k 






















1 

1 

1 

0 

0 

0' 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0^ 

1 

2 

2 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

1 

3 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

2 

4 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

3 

1 

5 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

3 

2 

6 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

4 

1 

7 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

4 

2 

8 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

5 

1 

9 ~ 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 
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16 145 OOOOOOOOOlOOOlOOOlOO' 

16 2 46 00000000010001000010 

16 3 47 00000000010001000001 

17 1 48 00000000010001001000 

18 1 49 *^0010000000010100000 

19 1 50 00010000000010100000 

20 1 51 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 

7.1. 2.1 Methodology, solutions and observations 

As mentioned in the previous section 7. 1.1.2 we do not use MWSC in the case of 

generalized grouping problem. Further, the proposed heuristic is applied in whole (i.e. 

all the steps) to the generalized grouping problems. However, in small size problem 

few steps may not get executed. The size of the problem is given by (plan x machine). 

With the Kusiak’s model we get the results as process families which is then 

transformed into corresponding part famihes. To reach at the optimal number of 

process families to be formed the model is run few times by trial and error (by varying 

the values of p). Machine assignment and refinement procedure of the proposed 

heuristic is also used after getting the solution fi-omthe formulation. 

With the proposed model (MAAM), we get the results as process families which may 

not be disjoint for bigger size problem (plans x machines > 50) and thus merging step 

(step2) of the proposed heuristic must be applied. After that machine assignment (step 

3) and refinement procedures (step 5)of the proposed heuristic are applied to assign 

machines to the famihes and to improve the solution, respectively. The final results are 

in terms of part famihes and their corresponding machines cells. 

Solutions of Example Problem 6 

Problem size : (11x4) and total number ofoperations(e') = 20 

Results by the three models are the same and shown below and in the Table 7.18. 

Total number of operation in the resulting solution (e) — 9. 

Number of bottleneck operations (Co) “ 0 and number of voids (ev ) ~ 1- 



Table 7.18: Solutions of various models using Lee’s distance Coefficient for example problem 6 


Model 

porl 

CPU 

time 

(sec) 

Process 

fiamilies 

Part 

femilies 

Machine 

cells 

Remarks 

Kusiak’s model ^MM) 

p = 2 

0.02 

_ (2,7) 
(5,9,11) 

(1,3) 

(2,4,5) 

(2,4) 

(1,3) 


Agarwal’s modified 
model 

p = 2 

0.09 

- do- 

- do- 

- do- 


Proposed model 
(MAAM) 

1 = 2 

0.03 

- do- 

- do- 

-do- 

no merging 
required 


Note : For small size problems (plans x machines <50), Kusiak’s model ma y take less 
time than the proposed model but the difference is just insignificant. 

Efficiency measures are : 


Efficiency measures M models 
grouping efficiency 95.00% 

grouping efficacy 90.00% 

global efficiency 100.00% 

group efficiency 100.00% 


Solutions of Example Problem 7 

Problem size : ( 1 5x5), total number of operations (e') = 34. 

The results obtained by all three models are the same and are diown below and in the 
Table 7.19. 

Total nmnber of operation in the resulting solution (e) =17. 

Number of bottleneck operation (intercell movement) (eo) = 1. 

Number of voids over all the part families (ev ) = 3. 


Table 7.19 : Solutions of various models using Lee’s distance coefficient for example problem 7 


Model 

porl 

CPU 

Process 

Part 

Machine 

Remarks 


time 

families 

families 

, cells 




(sec) 



(1,2,3) 


Kusiak’s model 

p = 2 

0.03 

(2,5,10, 

(1,2,5, 


(PMM) 



12,14) 

6,7) 

< 

I (4,5) 





(6,8) 

(3,4) 

• 

Agarwal’s 
modified model 

p = 2 

4.07 

- do- 

-do- 

3 - do - 


Proposed 
model (MAAM) 

1 = 2 

0.03 

- do- 

- do- 

: -do- 

merging 

required 



EflBciency measures are : 


Efficiency measures All models 


grouping efficiency 88.57% 

grouping efficacy 80.00% 

global efficiency 94. 1 1% 

group efficiency 88.88% 


Solutions of Example Problem 8 

Problem size : (5 1x20) with total number of operations (e') = 180 
The results obtained by all three models are shown in Table 7.20. 

Total number of operations resulting in all the solutions (e) = 67 


Table 7.20 : Solutions of various models using Lee’s distance coefEcient for example problem 8 


Model 

p or 1 

CPU 

Part 

Machine cells Co 

Cv Remarks 



time 

families 





(sec) 




Kusiak’s 

■ p = 5 

4.71 

(1,2,3,4,5) 

(1,6,9,12,20) 2 

17 no improvement by 

model 



(6,7.8, 

(2,5,7,16) 

refinement 

(PMM) 



. - . 9,10) 


procedure. 




(11,12,13) 

(3,8,11,19) 





(14,15, 

(10.14, 





16, 17) 

17,18) 



- - - ' 

- -- - 

(18,19,20) 

(4,13,15) 


Agarwal’s 

p = 5 

<= 

no integer 

- 

- 

modified 


50 

solution 



model 



possible 




Proposed 1 = 7 0.58 (1,2, 3,4,5) (1,6,9,12) 

model 

(MAAM) (6,7,8,9, (2,5,7,16, 20) 

10 ) 

(11,12,13) (3,8,11,19) 

(14,15, (10,14, 

16, 17) 17, 18) 

(18,19,20) (4,13,15) 

Efficiency measures are ; 

Efficiency measures 

grouping efficiency 
grouping efficacy 
global efficiency 
group efficiency 


Kusiak’s model and > Proposed 

AgarwaTs modified model : model 

95.25% 96.00% 

77.38% • 80.48% 

97.01% - 98.05% 

95.74% 97.87% 


1 15 best result is 

achievable. 

Cell size limit = 7 

large size problems 
are solved more 
efficiently. 



For studying the performance about the CPU time requirement, we take two more 
examples of OWRF type. They are given in the following tables. Observation on p or 1 
and CPU times have been given in Tables 7.21 and 7.22. 

Example Problem 9 (problem type OWRF) 

Fh^oblem size : (31x10). 

Number of parts Number of plans e' e 

10 31 131 42 


Table 7.21 : Observations on p or 1 and CPU-time for example problem 9 


Model 

p or 1 


CPU time (sec) 

Kusiak’s model (PMM) 

p = 3 


0.57 


p = 4 


0.15 

Agarwal’s modified model 

p = 3 


Disc quota reached 


p = 4 


( > 50 sec) 

Proposed model (MAAM) 

1= 5 


0.06 

Example Problem 10 (problem type OWRF) 

Problem size : (26x10) 

Number of parts Number of plans 

e' 

e 

15 26 


69 

47 

Table 7.22 ; Observations on p 

or 1 and CPU-time for example problem 10 

Model 

p or 1 


CPU time (sec) 

Kusiak’s model (PMM) 

p = 3 


0.13 

Agarwal’s modified model 

p = 3 


Disc quota reached 


p = 4 


( > 50 sec) 

Proposed model (MAAM) 

1 = 5 


0.05 



7.1.2.2 Graphical analysis of solutions 

Graph for number of plans, parts and e' vs CPU time for all the models is shown in Fig. 
7.4. 


figure 7.4 : Number of plans (or parts or e' ) vs CPU time 
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The conclusions here are the same as in case of simple grouping problems, that is, 
quadratic CPU-time variation for the proposed model (MAAM), and exponential 
CPU-time variation for Kusiak’s model. 


' Kusial 
Modifi 
— MAAN 


Graphs for various example problems vs four efiBciency measures for various models 
are shown in Figs. 7.5, 7.6, 7.7 and 7.8. 
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figure 7.5 iProblem number vs Grouping efficiency 
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FIGURE 7.7 ; Problem number vs Global efficiency 



FIGURE 7.8 : Problem number vs 

efficiency 


Group technology 

— □ Kusiak 
j - A - Modified 

■ • • X- - -MAAM 



For most of the examples taken (simple "as well as generalized grouping), aU of these 
measures are at least better for the proposed model (MAAM). 
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12 VALroATION OF PROPOSED SIMn^ARITY AND 
DISSIMILARITY COEFFICIENTS 

We test some example problems on simple grouping taken from the Uterature having 
process time and operation sequence inputs and the grouping results are obtained by 
applying their respective coefficients WPTSC and OSMUDC, in either Kusiak’s model 
(PMM) or the proposed model (MAAM). These results are then compared with their 
corresponding results from the literature. 

7.2.1 Observations on Weighted Part and Time Similarity Coefficient 
(WPTSC) 

Exa mple Problem 1 1 

Tlie problem by Wu and Salveudy [1993] is shown in Table 7.23. 

Table 7.23 : Process Time Part-Machine Incidence Matrix for Example Problem 11 

Machines 

1 2 3 4 5 6 7 

1 5 1 0 0 0 9 0 

P23007000 
a31380000 
r40006005 
t 50030700 
s 6 0 0 10 0 11 0 0 

7 7 2 0200 0 

Wu and Salvendy applied the network flow approach for rainimiTing the time of 
interceU movement between the machine cells. They obtained first the machine cells, 



and later the parts were allotted to the cells. The solutions obtained by them are given 
in Table 7.24. 


Table 7.24 : Wu and Salvendy’s solution for example problem 11 


Machine cells 

Part families 

(iZ) 

(2,4) 

OiQ 

(1,7) 

(2.3.5) 1 

(3,5,6) 

T otal number of intercell movements 

5 

Total time of intercell movements 

8 units 


Solution by using WPTSC (with weighing factors Ip = 5.0, h = 2.0.) in the proposed 
model (MAAM) is given in Table 7.25. As expected from MAAM, we obtain first the 
part families, and machines are then assigned to the families to form the machine ceDs. 
Table 7.25 ; Solution by using WPTSC with MAAM for example problem 1 1 


Part families 

Machine cells 

(1,3,7) 

(1,6,2) 

(4,2) 

(4,7) 

(5,6) 

(3,5) 

Total number of intercell movements 

2 

Total time of intercell movements 

13 units 


For this example, when we take Ip = 2.0, It = 5.0 we again find the same solution as 
given in the Table 7.25. 

Experiment 

Treating machines as parts and parts as machines, the result is obtained first in the form 
of machine cells to which parts are later allotted. Further, it is noted that applying 
WPTSC in this reversed case, we find exactly the same results as obtained by Wu and 
Salvendy (given in Table 7.24). Tims, we can say that WPTSC is an effective 
representation of similarity between a pair of plans, and also an effective measure to 
minimize the time of inter cell movements. 








Example Problem 12 

Ws problem has been solved by Seth [1995] using Mosier's process time related 
coefficients viz.. Additive Weighed Similarity Coefficient (AWSC) and Multiplicative 
Weighed Similarity Coefficient (MulWSC) developed by Hosier [1985], 


Table 7.26 : Input Data for Exanq)le Problem 12 
Machines 




1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

1 

1 

1 

1 

1 












0 

1 

2 

3 

4 

5 


1 

1) 

3 " 

0 

0 

0 

0 

0 

0 

0 

1 

3 

2 

0 

0 

0 


2 

0 

0 

4 

0 

2 

0 

0 

1 

0 

0 

0 

0 

3 

0 

2 

p 

3 

2 

0 

0 

0 

0 

2 

0 

0 

2 

0 

0 

0 

0 

2 

0 

a 

4 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

4 

0 

r 

5 

0 

0 

2 

0 

3 

0 

0 

3 

0 

0 

0 

0 

3 

0 

2 

t 

6 

1 

0 

0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

s 

, 7 

0 

4 

0 

0 

0 

0 

4 

0 

0 

4 

4 

4 

0 

0 

0 


8 

0 

0 

3 

0 

3 

0 

0 

3 

0 

0 

0 

0 

3 

0 

3 


9 

0 

0 

0 

3 

0 

1 

0 

0 

4 

0 

0 

0 

0 

3 

0 


10 

0 

3 

0 

0 

0 

0 

2 

0 

0 

4 

2 

1 

0 

0 

0 


Solutions by applying Hosier coefficients in Kusiak’s PMM model are given in the 
Table 7.27. 


Table 7.27 : Solution by using AWSC/MulWSC in PMM for Exan 5 )le Problem 12 

Machine cells 
(2,7,10,11,12) 

(3,5,8,13,15) 

(1,4,6,9,14) 


By using WPTSC (Ip = 5.0, h = 2.0) in the Kusiak’s model we obtained the same result 
as obtained by Seth (Table 7.27). Here also it is found that there is no effect of 
reversing the values of weighing factors Ip and It of part movement and time movement 
Thus, it can be concluded that MuWSC and TMDC are only different representation of 


Part families 

(2,5,8) 

13 , 4 , 6 , 9 ) 


the same kind of simUarity between the parts. However, it can be noted that this 
conclusion is problem specific. 

7.2.2 Observations on OSMUDC 

Examp lejn:oblejn 13 


The problem is shown in Table 7.28 with total number of operations (e) = 78 
Table 7.28 ; Input Data for Exanq)le Problem 13 


Number of parts 

20 

Number of machines 

20 

Part number 

1 

Sequence of machines required 
_ 12 1 9 18 20 

2 

" . 113 2 

3 

. 8 20 19 

4 

3 11 2 10 

5 

.4 15 6 7 

6 

"ll 14 16 17 5 

7 

7 5 16 17 

i" 

’ 15 _'l3 7 9 4 

9 

” 18 9 11 1 12 

10 

19 20 8 

11 

11 14 3 

12 

'9 18 5 12 1 

13 

6 7 15 17 

14 ^ 

8 10 12 

Is" 

‘l3 14 16 17 


15 7 6 19 

...... - 

9 1 12 

18 ~ 

8 19 20 10 

^19' 

3 2 11 5 

20 

18 10 1 12 


Seth [1995] obtained the solution by using SLBSC [Choobineh, 1988] in Kusiak’s 
model for part families formation and the machine cell formation in that order. The 


solution is shown in Table 7.29. 



Table 7.29 : Solution using SLBSC in Kusiak’s model 

Machine ceils Partfamflies 

.(.2,3An,13, 14,16, 17) (2,4,6,7,11,15,19) 

JQA1,?^L2J8} (12,1,17,9,20,14) 

i8,!9.20) . (3,10,18) 

(4,6,7,15) _ _ (5,8,13,16) 

Total number o f intercell movements = 10. 


EfiBciency measures are : 

Global efiSciency 87.18% 

Group ejQSciency 81.13% 

_^T. efficiency 75.86% 

Now, the results are obtained by using OSNxUDC (with inputs Wb = 1.0, Wo = 5.0 w; 
= 2.0, and UPm = number of operations under machine machine m) in Kusiak’s model 
(with p = 4) and applying steps 3 and 5 of the heuristic proposed are shown in Table 
7.30. 

Table 7.30 : Solution by using OSMUDC in Kusiak’s model 

Part famil ies _ ^ Machine cells 

(2,4 , 6,^iTT 5,19) " ^ (2, 3,5,1 1,14,16,17) ' 

17,20) “ (1,9,10,12,18) 

(3,To,18) (8,19,20) 

(4,6,7,15,13) 

Total number of intercell movements = 10. 

Efficiency measures are : 

Global efficiency 87.18% 

&oup efficiency 81.13% 

G.T.effiiciency 77.58% 

An improvement in G.T. efficiency is observed because machine 13 is kept in a 
different cell, however, that does not affect other efficiencies as the total number of 
intercell movements in both the cases remains the same. 


PHASE m 

7.3 TESTING OF PROPOSED IMPROVEMENTS IN KUSIAK’S 
MODEL (PMM) 

The testing is performed by solving some of the exanqrle problems of sec. 7. 1. 1. 1 using 
Lee’s distance coefiBcient in the improved models : MPGRP and NATGRP. The 
purpose of the in^rovement that at least two parts (plans of different parts) should 
always be included in a process family, is also observed. A comparison of CPU time 
taken by Kusiak’s model, EMPGRP model, and NATGRP model for given values of p 
and B (the maximum number of parts desired to be kept together in a femily) is also 
performed. We take example problems 1, 2, 3, 4 and 5 (of sinple grouping type) and 8 
(generalized grouping of OWRF type) for this purpose. 

7.3.1 Observation on Part Family Formation and CPU Time 

Observations on part family formation and CPU time are shown in Tables 7.31, 7.32, 
7.33, 7.35, 7.36 and 7.37. Table 7.34 shows the results of experimentation with 
varying the values of B on exan^rle problem 3. 

Observations on Example Problem 1 


Table 7.31 : Solutions of PMM, MPGRP and NATGRP models for example problem 1 


Model 

B p or 1 

CPU time (sec) 

Single part (plan) family 

Kusiak’s 

p = 4 

0.04 

Yes, 

model (PMM) 



part number 7 

IMPGRP ' 

B = 10 p = 4 

0.06 

No 

NATGRP 

B = 10 1=6 

0.08 

No 


Observations on Example Problem.? 

Table 7.32 : Solutions of PMM, IMPGRP and NATGRP models for example problem 2 


Model 

B 

p orl 

CPU time 

Single part (plan) 



^sec) 

family 

Kusiak’s model 


II 

0.02 

Yes, 

(PMM) 

IMPGRP 

' B = 5or l0 

p = 4 

0.01 

part number 2 
No 

NATGRP 

B = 5orl0 

1 = 4 

0.01 

No 


Observations on Example Problem 3 


Table 7.33 : Solutions of PMM. MPGRP and NATGRP models for example problem 3 


Models B porl CPU time Single part (plai 

- - - . . (sec) femily 

Kusiak’s model - p = 6 0.11 Yes, 

(PMM) _ part number 8 

IMPGRP B = 6 p = 6 0.11 No 

NATGRP B = 6 1=8 0.06 No 


Table 7.34 : CPlJ:time Experimentation with B in IMPGRP and NATGRP models 


CPU time (sec) 


B 

IMPGRP model 

NATGRP model 

4 

( > 50 sec) 

0.06 

5 

( > 50 sec) 

0.17 

6 

0.11 

0.06 

7 

0.12 

0.05 

F 


0.09 

9 

0.20 

0.09 

10 

*0.44 

0.24 

11 

0.26 

0.14 

12 

( > 50 sec) 

0.12 

13 


0.13 

14 


0.12 


7? 

0.10 

16~ 

7 ? 

0.13 

^ ”17 

77 

0.13 

18 

77 

0.06 

I? 

77 

0.13 

"""20 

77 

o.io 

21 

7 ? 

0.09 



77 

0.15 


We observe from the above table that for NATGRP model (ie. in case of natural 
grouping), the variation in CPU time with reqpect to values of B is negligible as 
compared to those with IMPGRP model (ie. in case of mqiosed grouping). Hence, any 



reasonable value of B can be used in the case of NATGRP model while for MPGRP 


model, a value of B needs to be explored for a minimum CPU time to result. This is so 
because of imposing the constrain of number of families in IMPGRP model 
Observations on Example 4 


Table 7.35 : Solutions of PMM, IMPGRP and NATGRP models for example problem 4 


Models B 

p orl 

CPU time 

Siugle part (plan) 



(sec) 

family 

Kusiak’s model 

p = 6 

0.25 

Yes, 




part number 6 

IMPGRP model B = 5,7 or 10 

p = 6 

( >30 sec) 

- 

NATGRP model B = 7' 

1=11 

7 0.07 

No 


Obseryations on Exam 


Table 7.36 : Solutions of PMM, IMPGRP and NATGRP models for example problem 5 


Models 

B 

porl 

CPU time (sec) 

Single part (plan) 





family 

Kusiak’s model 

- 

,,P.=,7 . 

9.41 

No 

IMPGRP model _ 

"¥=’16^ 

P,= 7 

( > 50 sec ) 

- 

NATGPyPmodd^ 

71^= i63'„ 

1= 18 

0.98 

No 


Observations on Example ^ 


Table 7.37 : Solutions of PMM, IMPGRP and NATGRP models for example problem 8 


Models 

B 

p orl 

CPU time 
(sec) 

Single part (plan) 
femily 

Kusiak’s model 


p = 5 

5.44 

No 

(PMM) 


P = 7 

18.93* 

No 


“ 10 

p = 5 

( > 50 sec) 




, P = 7 . 

2.60 

No 

NATGRP modef^ 


1=8 

1.17 

No 


' ' 10 

1 = 8 

1.81 

No 


Larger value of ‘p’ in Kusiak’s model not necessarily always results in lesser 
CPU time as can be observed in the above table. 



CPU time (sec) 


7.3.2 Graphical Analysis of Solutions 

Graph of total number of plans and e' vs CPU time is shown in Fig. 7.9. 


10 


0 - 
8 - 
7 - 
6 - 
5 - 
4 - 
3 - 
2 - 
1 - 


0 th 
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It is found that with NATGRP model, CPU time varies quadraticaJly with the total 
number of plans while that with Kusiak’s model it varies e?q)onentially. However, 
IMPGRP is not able to solve the problem of large number of plans within permissible 
time. These observations validate their conoparisons already given in the Table 4.2 on 
the basis of the total number of enumerations. 

Some interesting findings about the similarity between the pair of models are worth to 
be noted. It is seen that NATGRP and the proposed models are somewiat similar with 
respect to solutions CPU times and the final solutions as they both performs natural 
grouping. The CPU times for the case of NATGRP model are, however, little higher 



than that for the case of proposed MAAM model. This is due to larger number of total 
constraints (2q + K) in case of NATGRP model than that in the case of the proposed 
model (2K + q) [see Table 4.4], while Agarwal’s modified model and IMPGRP model 
are similar in the way of grouping i.e. inq)osed grouping and very hi^ CPU time, their 
differences are the violation of GT concept by AgarwaTs modified model and the value 
of B. This value of B is the maximum niunber of parts which are desired to be kept 
together in a family for the IMPGRP model and (K-p+1) for Agarwal’s modified 
model. 



8 


CONCLUSIONS AND SCOPE OF FUTURE 
RESEARCH 

7.1 CONCLUSIONS 

In the present work, a new integer programming model (IMAAM) has been presented 
for clustering parts to form part famihes in the generalized grouping environment 
where each part can have alternative process plans and their operations can be 
performed on alternative machines. The proposed model (MAAM) is shown to be 
equivalent to unit capacitated circulation network flow (UCCNF) model with the 
objective of minimizing the cost of flow. The resulting solution is in the form of a main 
loop and several sub-loops which are interpreted as process femilies. The model does 
not require any knowledge about the number of femilies to be formed a priori vshich is 
one of the requirements in the case of p-median model (PMM) of Kusiak[1987]. The 
model ensures that at least two parts will be contained in a femily while Kusiak’s 
model does not. Singleton family i.e. family of single part may get formed by Kusiak’s 
model. The model produces the number of famihes by itself However, as the families 
formed are large in number and also not disjoint of each other for large problem, so, a 
simple heuristic procedure is proposed to reach to the optimal solution. In feet, part of 
the heuristic is common to the Kusiak’s model also. The additional steps of the 
heuristic which is required for the resulting solution fi-om the proposed model 



(MAAM) is of <9(1 ), where 1 is the number of loops resulting in the solution. The 
proposed model itself has a low order of complexity that is of quadratic order, in terms 
of total number of feasible enumerations. Hence the overall conq)lexity is also of 
quadratic order only, while for Kusiak’s model it is of the exponential order because of 
the prespecified number of part-families. Agarwal’s modified model has also been 
included for the purpose of comparing the results with the proposed model. It has been 
observed that the proposed model (MAAM) is computationally extremely efficient as 
compared to both Kusiak’s model and Agarwal’s modified model. The proposed 
model (MAAM) gives the optimal solution in one go while Kusiak’s models need to be 
solved few times. Agarwal’s modified model is not able to solve even medium size 
problem in the permissible time. It has been found that the heuristic procedure when 
applied to the solution resulting firom the proposed model (MAAM) results into the 
best achievable solution with respect to the minimization of intercell movements. In 
case of Kusiak’s model, application of the refinement procedure of the heuristic may 
result into the best achievable solution depending upon the preq)ecified number of 
families to be formed. If it is what comes as a result of applying the proposed model 
(MAAM) and the heuristic, then the resulting solution of the Kusiak’s is the best 
achievable one. 

Various efficiency measures are calculated for the solutions of these models. It has 
been found that for the proposed model (MAAM), these measures responded at least 
better than other models, for most of the problem taken fi:om the literature. 
Improvements in Kusiak’s model have also been earned out. These are . Inqrosed 
grouping model (IMPGRP) and natural grouping model (NATGRP). They ensure that 
at least two part is contained in a part femily t’^duch is otherwise violated by Kusiak s 



model. However, we require ‘B’, the maxiinum number of parts making a part femily, 
as an input in these improved models. This number has a major effect in MPGRP 
model foi the solution to result in the least conq)utational time, but for NATGRP 
model it is almost ineffective. The NATGRP model is found to have similarity with the 
proposed model (MAAM). This has also low order of complexity (see Table 4.2) and 
the resulting solution need to go through the proposed heuristic. MPGRP also need to 
go through some steps of the heuristic proposed. 

Some measures of similarity and dissimilarity have been developed to arrest the 
similarity between the plans having inputs other than 0-1 nature, like, processing time 
of operations and their sequence. Some examples of the literature have been tested and 
compared with their existing solutions. Thus the validation of the developed 
coefficients is carried out. 

7.2 SCOPE OF FUTURE RESEARCH 

The proposed model (MAAM) can be tested with larger si 2 e problem including say 
500 machines and thousands of parts. The problem can be created by random number 
generation programming. Generalized grouping with operation sequence can be 
viewed in such a way that the operation sequence itself can be included in the network 
configuration. The proposed model can be extended considering the required number 
of part families to be known in advance. A clue in this direction can be given as 
extending the dimension of variable Xy of 2 to 3 i.e., Xyk, where k is the part number. A 
new scope can be the consideration of capacity constraints of machines in either some 
type of measures or in the formulation itself The clustering problem under random 
product demand is a new avenue for research. As the GT problem is NP-hard, so some 



robust heuristic approaches like Genetic Algorithm, simulated annealing method, tabu 
search etc. can well be examined. Other new tools like fuzzy logic, neural nets can also 
be explored for their application in group technology problems. 
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appendix a 

Explanation of Constarints for Example Problem 

Tlie constraints for the example problem given in Table 4. 1 for the case of the 
proposed model (MAAM) are explained as below. 

Con straints of eq. (4.4) 

Xi 4 + Xi 5 Xif, Xi 7 + Xi 8 + Xi9 + Xno + Xni + X24 + X25 + X26 + X27 + X28 + X29 + X210 + 
X21 1 + X34 + X35 + X3f. + X37 + X38 + X39 + X310 + Xsn = 1 

X4I + X42 + X43 X4(, + X47 + X 48 + X49 + X 410 + X4n + X51 + X52 + X53 + X56 + X57 + X58 + 

X59 + X.MO + X 5 n “ 1 

x<,i X(,2 X<,.1 ‘t X<,4 ^ X <,5 + Xfis + Xfit) + X^io + Xen + X71 + X72 + X73 + X74 + X75 + X78 + 

X79 t X 7 U) * X71, - 1 

Xki t Xk2 ♦ Xh 3 f X 84 * Xk? 4 ' Xh(, •+■ X 87 + Xsio + Xsn + X91 + X92 + X93 + X94 + X95 + X95 + 
Xt )7 XiJlO X <)11 1 

XlOl t Xi ()2 < Xi ()3 f Xi().i -t Xi ()5 •+ Xi 06 + X107 + X ]08 + X109 + Xnoi + X1102 + Xn 03 + Xno4 + 

X|U )5 * Xiiw. * Xiio? < XilOK + Xn ()9 = 1 

Const rainl s of cq. (4 .5) 

X41 ^ X51 f x<,i ) X71 ♦ Xki ■+ X91 + Xioi + Xiioi + X42 + X52 + X62 + X72 + X82 + X92 + X102 + 

XU02 ♦ X43 t X53 + Xa ’+ X73 + x «3 + X93 + Xi 03 + Xn 03 = 1 

Xi4 ■< X24 t X<4 -t X<4 + X74 ■+' Xk 4 + X94 + Xi04 + X1104 + Xi4 + X25 + X35 + X65 + X75 + Xss + 

X95 ^ Xi ()5 ^ XiiOS 1 

Xl(, X26 t X36 t X 4 <. t X 56 + Xsfi + X% + X106 + Xji 06 + Xn + X27 + X37 + X47 + X57 + X 87 + 
Xy? Xi()7 "t Xij()7 ■ 1 

XlK ■+ X2K + X 3 K ^ X 4 K X 58 + XfiK + X78 + X1O8 + Xn08 + X19 + X29 + X39 + X49 + X59 + X 69 + 

X79 XuK) + Xu(W ^ 1 

Xiio + X210 + XiKi 1 X410 + X510 + x<iio + X710 + Xsio + X910 + Xni + X2n + Xsn + X 4 n + Xsn 
■+'X^,u+ X 7 U + Xki I + Xyn = 1 

Con.straint.s of eq. (4.6) 

Xi 4 t X 4 J t Xi5 + Xsi + Xi 6 + X 6 I + Xi 7 + Kn + Xi 8 + X 8 l + X19 + X91 + Xuo + Xioi + Xm + 

XIKIJ 0 or 2 

X24 t X42 -t X25 i X52 + X26 + X4,2 + X27 + X72 + X28 + X82 + X29 + X92 + X21O + X102 + X2n + 

Xuo2 0 or 2 

X34 + X43 f X35 + X53 + X36 + Xe3 + X37 + X73 + X38 + X 83 + X39 + X93 + X310 + Xi 03 + X311 + 

Xno3 0 or 2 

X41 * Xj4 ■+ X42 + X24 + X43 + X34 + X46 + X64 + X47 + X74 + X48 + X 84 + X49 + X94 + X 410 + 

Xi()4 + X4ii * xno4 0 or 2 

Xf.l Xk, ■( 50,2 < X2<, + x «3 + X 36 + X64 + X46 + XfiS + Xs6 + X68 + X86 + X 69 + X96 + XfilO 

Xw, + X4,ir^ xihk. 0 or 2 4. -i- 

^71 + )C}7 4 %72 4 K27 ^73 ^ ^37 ^74 + X47 ^57 7C78 + Xg? ^79 X97 X710 

Xi()7 + X71J 4 xno7 =” 0 or 2 
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txx:~ 0 :: ; ^ ^ "» ^ ^ ’‘» ^ ^ ^ ^ - x... - 

Xio, + X „0 + X,„2 + X2,o + X,„3 + X3.0 + X,„ + X 4,0 + X,„; + Xj.o + X,« + X„„ + X,„, + X,,„ 

+ X108 + Xkio + X)o') + X910 = 0 or 2 0 710 

X,,., + X.„ + X ..„3 + X3„ + X„« + X3„ + X,.„. + X 4 „ + X „.3 + X 3 ,. + X„<« + X.„ + X„o, 
+ X 7 U + XjiOK + X«u + Xu 09 + X 9 U = 0 012 

Thus, total number of constramts = 2K + q = 2x5 + ll = 21 
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